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1. Introduction 

To high energy theorists, a chiral gauge theory is not an unfamiliar concept. A 
Dirac spinor consists of two Weyl spinors of opposite chiralities: the left- and right- 
handed ones. A chiral gauge theory is a theory of Weyl fermions, not necessarily 
forming Dirac multiplets, but living in some complex representation of a Lie group 
and minimally coupled to the corresponding gauge field. We will adopt a broad 
definition of a "chiral gauge theory" — a theory, where gauge- invariant mass terms 
for all fermions are not allowed — which includes some interesting cases with pseudo- 
real representations, where fermion mass terms vanish due to statistics. 

According to our experimentally verified knowledge, nature is described by just 
such a theory. It might sound a bit mysterious that we do not yet have a method of 
approximating an arbitrary chiral gauge theory by latticizing and then simulating 
it on a computer — even in principle. If one wishes to use the lattice to gain some 
insight into its non-perturbative properties, we just do not know how to do it. 

Apart from interest in physics of the Standard Model — which, at low energies, 
is a weakly-coupled spontaneously broken chiral gauge theory that does not obvi- 
ously call for a lattice study — interest in strong chiral gauge dynamics has both 
intensified and abated during the past few decades. From the overview in the next 
Section, it should be clear that while there exist potential applications of strong 
chiral gauge dynamics to particle physics, at the moment it appears difficult to 
identify "the" chiral theory most relevant to particle physics model-building (apart 
from the weakly- coupled Standard Model, of course). Thus, the problem of a lattice 
formulation of chiral gauge theories is currently largely of theoretical interest. This 
may or may not change after the LHC data is understood. Regardless, we find the 
problem sufficiently intriguing to devote some effort to its study. 

a Leading experts' published opinions on this matter vary, e.g.: "If a solution to putting chiral gauge 
theories on the lattice proves to be a complicated and not especially enlightening enterprise, then it 
probably is not worth the effort (unless the LHC finds evidence for a strongly coupled chiral gauge 
theory!)." (D.B. Kaplan^, or: "Without a proper lattice formulation of a chiral gauge theory, it 
is unclear whether such models make any sense as fundamental field theories." (M. Creutz^l. 



Chiral lattice gauge theories via mirror- fermion decoupling: a mission (im)possible? 3 



1.1. Chiral gauge theories through the ages 

Strongly coupled chiral gauge theories were hrst studied in the 1980's for their pos- 
sible relevance to quark and lepton compositeness. To illustrate both the connection 
to such ideas and the quite different physics that can occur in chiral gauge theories 
(compared to vectorlike, or "QCD-like" theories), consider the classic example^ of 
a four-dimensional strongly-coupled chiral gauge theory. This is an SU(5) gauge 
theory with a 5* and a 10 Weyl-fermion representations, both taken left-handed. 
The theory has one anomaly- free chiral U(i) global symmetry^] under which the 
5* fermion has charge —3 and the 10 has charge 1. It is a simple matter to check 
that the 't Hooft anomaly matchin gj 4 l 5 l 6 l conditions for the U(l) 3 and U{1) anoma- 
lies are saturated by a single massless composite Weyl fermion, a 5*-5*-10 bound 
state. Thus, in the absence of gauge-invariant order parameters to break the chi- 
ral symmetry, the likely infrared (IR) dynamics is that of a free single massless 
gauge-singlet composite fermion (a "tumbling" scenario based on "most-attractive- 
channel" (MAC) reasoning is consistent with such an IR spectrurrFJ). Since the 
anomaly-free chiral symmetry is unbroken, this chiral gauge theory exhibits "con- 
finement without chiral symmetry breaking" — clearly, an IR dynamics different from 
that of QCD-like theories (at least from vectorlike theories without supersymmetry 
and the associated massless elementary scalars: it turns out that supersymmetric 
QCD, a vectorlike theory, has phases that exhibit confinement without chiral sym- 
metry breaking, see^for an overview and references). 

The appearance of massless composite fermions in many chiral gauge theories, 
as suggested by anomaly matching and the (non rigorous) MAC arguments gave rise 
to many attempts to build "preon" models for elementary particle physics; see^ for 
references. This direction is currently not being actively pursued (at least not via 
direct studies of four-dimensional strong gauge dynamics, see below), not because 
compositeness is ruled oul]^] but because constructing realistic models turned out 
to be difficult and because the theoretical tools available, at least for nonsupersym- 
metric models, are limited to the already mentioned 't Hooft consistency conditions 
and MAC. However, the idea that some of the particles we see may be composite is 
still alive. More recently, after the advent of the AdS/CFT correspondence^, mod- 
els of electroweak-scale physics have used quarks, leptons, and gauge fields "living" 
at various locations in slices of AdSs. These phenomenological models are thought 
to be holographically dual to strongly-coupled gauge theories, where some of the 
Standard Model particles arise as (partially) composite objects^. However, it is 
currently not known what are the four-dimensional gauge theories whose dual de- 
scription is given by these phenomenological warped models, much less whether 
they are chiral or not. 

Another surge of interest in chiral gauge dynamics occurred in the supersymmct- 

b The quadratic Casimir is 1/2 for the fundamental and 3/2 for the antisymmetric tensor of SU (5). 
Experimental bounds are stronger than in the 1980's, constraining the scale of compositenessi^l 
to be in the multi-TeV region. 
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ric context in the mid-1990's. It was due to the then-held belief that only chiral gauge 
theories can be used to dynamically break supersymmetry. Recall that if supersym- 
metry is relevant to particle physics, it has to be broken — hence, the mechanism of 
supersymmetry breaking is of interest to superparticle phenomenology. A number of 
chiral supersymmetry-breaking theories were studied, see^ for a review, using the 
newly available theoretical tools of holomorphy and duality. The interest in chiral 
supersymmetric theories has somewhat subsided after the recent discovery, reviewed 
irP, that vectorlike supersymmetric theories can have metastable supersymmetry- 
breaking vacua, sufficiently long-lived for particle physics model-building purposes. 

Electroweak symmetry breaking and fermion mass generation is another area 
where strong chiral gauge dynamics has been invoked- 15 ! 16 !, for example usin^^ tum- 
bling (via MAC) in chiral gauge theories to construct extended technicolor models. 
A more recent rather speculative proposal 1 ^ is to break electroweak symmetry in a 
chiral theory of massless fermions (a fourth generation) with both electric and mag- 
netic U(l) charges. Consistent theories of "mutually nonlocal" massless monopolcs 
and dyons have been argued to exist in J\f=2 supersymmetry®, but it is not known 
whether nonsupersymmetric versions exist; addressing this on the lattice is chal- 
lenging, as such theories are inherently strongly coupled. It is difficult to assess 
the viability and study detailed features of the proposed strong chiral electroweak- 
breaking models, without credible tools to study their dynamic^] 

Finally, on the theoretical side, there has been little progress (in the nonsuper- 
symmetric context) in the nonperturbative understanding of chiral gauge dynamics 
since the 1980's^. An exception is the recent worUHZlMHll on studying the nonper- 
turbative dynamics of chiral gauge theories compactified on i? 3 x S . At sufficiently 
small radius of S 1 , semiclassical methods can be used to study the generation of 
mass gap in the gauge sector (i.e., the string tensions) and the breaking of discrete 
or abelian chiral symmetries. These studies show, once again, that nonsupersym- 
metric chiral gauge dynamics is different from that of QCD-like theories. However, 
while instructive, the validity of the semiclassical approach does not extend to the 
physical case of decompactified R 4 . 

1.2. Where does the lattice stand? 

Ideally, lattice studies should provide both a definition and a means to study the 
dynamics of any strongly-coupled gauge theory. In vectorlike nonsupersymmetric 
theories, significant theoretical and practical progress has been made in the past 15 
years. In chiral theories, there are both conceptual and technical obstacles. These, 
along with some of the ongoing attempts to overcome them, arc the subject of this 
review. 

We note that there has been some recent progress, reviewed irP3 towards a 

d Dcspitc the new poorly-understood strong chiral interactions at the TeV-scale, characteristid221i£l 
signatures of these models could be found at the LHC or their absence used to rule them out. 
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lattice definition and even numerical studies of some of the theories of possible 
phenomenological interest mentioned in Section most notably the ones with 
supersymmetry. Practical progress is, however, limited mostly to vectorlike super- 
symmetric theories with extended supersymmetry in lower dimensions; see^ for an 
exception. The chiral supersymmetric case, in both two (2d) and four (4d) dimen- 
sions is currently out of reach. As already stated above, for nonsupersymmetric 
chiral theories, the lack of theoretical tools to study such theories makes it difficult 
to assess the various models. 

The long-standing issue of chiral lattice gauge theories has been previously re- 
viewed ir i 27 ! 28 ! 29 !, where different approaches and their status are discussed in more 
detail than we can go into here. We also recommend the already quoted lectures^ on 
exact chiral symmetry in vectorlike lattice theories; we will not review this subject 
in great detail and will largely assume the reader's familiarity with it. 

1.3. Outline 

(1) In Section [2j we describe the problems that must be overcome in order to define 
chiral symmetry and chiral gauge theories on the lattice. 

• In Section |2.1[ we recall why the partition function of a chiral gauge the- 
ory is only defined up to a phase, a property not usually relevant in the 
continuum, but important for the lattice definition. 

• In Section [2~2} we briefly discuss the complexity of chiral fermion determi- 
nants (one of the main differences between chiral and vectorlike theories 
from the point of view of the Euclidean path integral). 

• In Sections |2.3| and |2.4| we recall the fermion doubling problem, its solu- 
tion in terms of "Ginsparg- Wilson fermions," and some of their pertinent 
features. 

• In Section |2.5| we list the various problems arising when using Ginsparg- 
Wilson fermions to define chiral lattice gauge theories. 

(2) In Section [3) we give a non-technical discussion of the "mirror-decoupling" 
approach to chiral lattice gauge theories. 

• In Section [3. 1[ the main idea is described using continuum notation, along 
with a discussion of some motivational evidence. 

• In Section |3.2[ we explain the nature of the strong-coupling symmetric 
phases in lattice theories with strong multi-fermion or Yukawa interactions 
and illustrate the desired features of the mirror-decoupling dynamics with 
the help of a simple solvable toy model. 

• In Section |3.3[ we use 't Hooft anomaly matching arguments to explain 
why breaking all mirror chiral global symmetries, anomalous or not, is a 
necessary condition for mirror-fermion decoupling. 

• In Section |3.4| we compare our approach to "pre-Ginsparg- Wilson" at- 
tempts of mirror /doubler decoupling. We emphasize that the crucial dif- 
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ference lies in the absence of exact lattice chiral symmetries in earlier 
implementations. 

• In Section 3.5 devoted to the "big picture," we end the more descriptive 
part of the review by giving a list of questions about the mirror-decoupling 
dynamics that need to be answered and a summary the (known and un- 
known) answers. 

(3) In Section |4j we study the many technical issues that arise in the Ginsparg- 
Wilson-fermion implementation of "mirror-decoupling" approach. 

• In Section |4.1[ we study general properties of chiral partition functions. 
Many theoretical results, more generally valid, will be illustrated on the 
example of a simple model, as discussed below. The tools developed allow 
the questions listed in Section 3J3 to be addressed, here and in future work. 

— In Section |4.1.1| we introduce the model that we have studied ex- 
tensively, the "1-0" 2d model. Here, the "mirror" fcrmions are in an 
anomalous representation, so we do not expect to decouple them with- 
out breaking the gauge symmetry. However, this model provides a 
useful playground to study the dynamics of strong Yukawa or mul- 
tifermion interactions with Ginsparg- Wilson fermions via inexpensive 
numerical simulations and allows addressing some of the questions of 
Section 13.51 

— In Section |4.1.2| we explain how to separate the action and measure 
of a vectorlike theory into "light" and "mirror" parts and discuss the 
"splitting theorem" about the dependence of chiral partition functions 
with general interactions on the gauge background. This is an impor- 
tant analytical result, without which the studies discussed here would 
have not been possible. 

— In Section [4. 1.3| we use the split partition function to study the phase 
structure of the mirror dynamics of the "1-0" model at strong Yukawa 
coupling. 

— In Sections 4.1.4 and 4.1.5 we study the mirror polarization operator 
in vanishing gauge background. Its real part is a universal probe of the 
mirror-fermion spectrum, while its imaginary part contains (in 2d) the 
mirror-fermion anomaly. We list some general exact (i.e., independent 
of the nature and strength of the mirror interactions) properties of 
the mirror polarization operator, useful for its analytic and numeric 
studies. 

• In Section |4.2| we explain how to formulate and study 't Hooft anomaly 
matching on the lattice, applied to the case of (strong or weak) mirror 
dynamics. 

• In Section |4~3} we discuss the results of the numerical study of the realiza- 
tion of anomaly matching in the mirror sector of the "1-0" model, both in 
the strongly-coupled symmetric and weakly-coupled "broken" phase. We 
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show that in the strong symmetric phase anomaly matching is satisfied 
with the minimal number of massless charged fermions (and stress the 
role of the Majorana coupling in ensuring minimality). 
• In Section |4.4| we rewrite the partition function in a different basis, mo- 
tivated by the relation of Ginsparg- Wilson to domain-wall fermions. This 
representation is useful to further understand aspects of the dynamics: we 
explain the appearance of extra massless vectorlike fermions, found in Sec- 
tion 4.3 when one of the mirror-theory couplings (the Majorana coupling) 
vanishes. We also comment on the relation of our formulation to other 
proposals to use strong coupling in the domain wall set-up to decouple 
mirror fermions. 

(4) In Section [5] we conclude with a brief discussion of the outlook for the future. 

2. Difficulties faced by chiral gauge theories on the lattice 
2.1. An ambiguous phase 

Before we explain why defining lattice chiral gauge theories is so hard, we should 
first mention that the partition function of any chiral theory is, rigorously speaking, 
not defined. This is because an operator that maps between two independent vector 
spaces does not have a natural definition of its determinant. Take for example the 
kinetic term of a pair of left-handed Weyl fermions: 

S = J d 4 x^ L Dijj L (1) 

where D stands for the usual Weyl operator and ipL an d ipL are two independent 
left-handed Weyl fermions (in Euclidean space, the left-handed spinors and their 
conjugates transform in different representations of the Euclidean "Lorentz" group 
50(4) ~ SU(2) x SU(2), not related by complex conjugation as in Minkowski 
space). Attempting to define the Grassmann path integral: 

Z = J D^ L Di; L e s , (2) 

one imagines choosing a set of orthonormal basis vectors, say Ui and u,, and ex- 
panding in them the fermion fields w ith Grassmann coefficients (ip L is expanded in 



u\ and ipL — in v%, see Section 4.1.21). The subscript i runs from 1 to the dimension 



of the fermion phase space, say d, which we for simplicity assume to be the same for 
ip L and ipL. Certainly d is infinite for theories defined in the continuum spacetime, 
and this discussion is rather formal at the moment. The partition function ^ is 
evaluated by just defining Z = det(u\Dvj), where u\Dvj is a dx d matrix (of indices 
i and j), whose determinant is always defined. Notice, however, that the bases {ui} 
and {vi} are not unique. We can equally well choose, instead of, e.g., v il a different 
set of vectors v[ = UijVj, where Uij is a unitary matrix. If we compute Z as above 

we find the result differs by a factor of det . This 
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factor is always a pure complex phase since [Ty is unitary. The argument presented 
here is equivalent to the trivial fact that the identification between two vector spaces 
of the same dimension is not unique. 

Such a phase ambiguity of the partition function is usually not a severe problem 
since it is sufficient to just choose and then stick to a particular set of basis vectors. 
The ambiguous phase of Z — provided it can be chosen independently of the gauge 
field, as in perturbative continuum regularizations — is always divided out in the 
vacuum expectation values of any operator, which are the physical observables. In 
what follows, we will see that in attempts to define chiral lattice gauge theories this 
ambiguous phase can not be taken to be gauge-field independent. Its determination 
constitutes the so-called "measure problem," that will play an important role in our 
further discussion. 

2.2. The sign problem in chiral gauge theories 

Even if the "measure problem" of Section |2.1| is solved, an important issue that 
remains is the sign problem of chiral partition functions. The sign problem is likely 
to hinder an immediate application of any successful lattice definition of chiral 
gauge theories to studies of strong gauge dynamics via Monte Carlo methods (as it 
generically occurs in 4d chiral gauge theories). However, as we discuss at the end 
of this Section, there are some exceptions — the continuum fermion determinant is 
real in some special cases in 4d as well as in 2d anomaly-free chiral gauge theories. 

From the point of view of the Euclidean path integral, generic chiral and vector- 
like theories differ in that in the vectorlike case the fermion effective action is real, 
while in the generic chiral case it has an imaginary part. The real part of the fermion 
effective action in a chiral gauge theory is equal to one-half that of the vector theory 
obtained by adding massless mirror fermions to the chiral theory. The appearance 
of an imaginary part in chiral theories can be explained in the continuum and we 
will briefly do so in this Section. In the Euclidean path integral language, the imagi- 
nary part[^]is the one that must be responsible, along with the distinct effect^HHl of 
nontrivial topological sectors, for the quite different dynamics of strongly-coupled 
vectorlike and chiral gauge theories (as seen in our discussion of the SU(5) chiral 
model in Section |1.1[ ) . 

The reason for the sign problem is that chiral gauge theories — such as the Stan- 
dard Model — are left-right asymmetric, i.e. violate parity. In Euclidean space, the 
parity-violating part of the effective action for fermions is purely imaginary. One 
can see this in perturbation theory by recalling that the fermion determinant can 
be represented by a sum of one-loop fermion graphs with an arbitrary number of 
gauge field insertions. The real, parity-even, part of the effective action sums over 
loops of vectorlike fermions with the quantum numbers of the left- and right-handed 

e To avoid confusion, we note that after integration over the gauge fields, the partition function is 
real — the imaginary contributions to the partition function of field configurations that are parity 
images of each other cancel (while the real parts are equal) leaving only a sign problerrPHl. 
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fermions of the chiral gauge theoryj^] 

1 00 1 

Re log Z[A]± = -- £ — Tr (D^iA^) 2n + c.t., (3) 
n=l 

where Do is the free Dirac-fermion massless propagator (we give the effective ac- 
tion for a single Weyl fermion; the expression ^ has to be summed over all Weyl 
fermions of the theory). The imaginary part of the effective action for a left- or 
right-handed Weyl fermion, on the other hand, is given by the parity-odd expres- 
sion: 



1 1 

1 Im log Z[A} ± = T - Tr 75 (D^iA^) 2n+1 + c.t. . (4) 

n=0 



Both <[3j) and @ are valid up to the addition of local counterterms, denoted by 
"c.t.," which ensure that the effective action is finite and gauge invariant (if gauge 
anomalies cancel). 

In 4d chiral gauge theories, the n > 1 terms in the imaginary part of the effec- 
tive action are finite and gauge invariant — and must contain important information 
about the chiral gauge theory dynamics, as we stipulated above. Using Pauli-Villars 
or zeta- function regularization, one can show^that the difference between the imag- 
inary part of the fermion effective action Q computed for two different perturbative 
gauge backgrounds, A and A Q , in the anomaly- free case, can be written as: 

Im \ogZ[A] - Im \ogZ[A ] = vr^O] . (5) 

Here, 77 [0] is the "77-invariant," or spectral asymmetry, of a five-dimensional her- 
mitean Dirac operator H constructed from the the four-dimensional vectorlike Dirac 
operator (A^ — A^T a is hermitian): 

H = 75 i^ + iD[A(t)], D[A] =D - ij^ , (6) 

where A(t) smoothly interpolates between A{— 00) = A and A(oo) = A. The 77- 
invariant 77 [0] appearing in ([5| is the analytic continuation to s — of 77 [s], denned 
for large Re(s) > as: 

n\s] = Tr s g . (7) 

The definition has to be slightly modified when zero modes of H are present.^ 

The 77-invariant representation of the imaginary part of the chiral determinant 
is useful to study anomaly cancellation, as it can be shown that the change of the 

f To obtain the formal expressions (J3J and one expands logZ[A]±=logdetD± (with D±=Dq — 
ijfiA/j, 1 2 5 1 = 7^9^) in a perturbative series and separates the real and imaginary parts; we 
use hermitean Euclidean gamma matrices. We note also that similar expansions for the real and 
imaginary parts can be obtained using the lattice definition of the determinant via the Neuberger- 
Dirac operator; naturally, these coincide with eqns. J3l, Q if the naive continuum limit is taken 
first. 
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effective action ^ under gauge transformations of A is a topological invariantl 32 ! 33 ! 
On the other hand, the imaginary part of the action itself is not a topological term 
and relating it to the 77-invariant, computed with some fiducial background field, 
e.g., Aq in ([5]), sheds no light on dynamical aspects of the physics. 

Another use of equations ^ and Q is to quickly argue that the imaginary part 
of logZ[A] vanishes for real representations. This follows from the fact that for real 
representations, there exists a gauge group element a such that (iA^)* = a~ 1 iA ll a. 
If \& is an eigenfunction of H, eqn. (6|, with eigenvalue A, then 1 ff' = B~ 1 a~ li ff* 
with B~ 1 j /M B = 7* has eigenvalue^ 



-A. Consequently rj [s], eqn. (|7j), vanishes for 
any s, and so does the spectral asymmetry rj [0] and the imaginary part of the 
fermion effective action. 

The result of the previous paragraph is of interest for some chiral gauge theo- 
ries. There exists a 4d chiral theory that has a real positive determinant — the chiral 
SU{2) gauge theory with a single Weyl fermion in the three-index symmetric repre- 
sentation. Despite the pseudo-reality of the j = 3/2 representation, a fermion mass 
term vanishes by Fermi statistics and the theory is usually referred to as chiralj^] 

In 2d anomaly-free chiral models, the terms with n = in the imaginary part 
Q cancel, while terms with more than two powers of the gauge field in the Weyl de- 
terminant (the sum of ([3J and Q) can be argued to vanish^ due to their regulator- 
independence, gauge invariance, and chirality. Thus, at least in the continuum, 2d 
anomaly-free chiral gauge theories can be defined to have a real determinant. That 
this is so becomes explicit upon examining their solution via bosonizatiorE^EH 

We expect that the positivity of the fermion determinant should be reproduced 
in the continuum limit by a lattice definition, hence (assuming that this implies 
a less severe sign problem on the lattice) a successful proposal for defining chiral 
lattice gauge theories stands a chance to be tested via Monte Carlo simulations with 
dynamical gauge fields only in 2d or in some special 4d chiral gauge theories, like 
the j = 3/2 SU(2) theory mentioned above — at least in the foreseeable future. 



2.3. The fermion doubling problem 

As a most simple-minded attempt to define lattice theories of fermions, one may just 
discretize the Dirac or Weyl operator as one would normally do for the Laplacian. 
Not surprisingly, this does not work, or we would not have to write this review. The 
problem is that when one takes the continuum limit of a lattice theory with a naive 
discretization of the Dirac or Weyl operator, one finds, most intriguingly, that every 
Weyl fermion introduced in the action is accompanied by an unexpected massless 
doubler of opposite chirality. 

s This follows from aBHB~ 1 a~ 1 = —H*, which is easy to check. 

h The supcrsymmetric version of this theory has been used as a model for dynamical supersymmetry 
breaking, but recent argument£2E3 suggest that it is conformal in the IR; the nonsupersymmetric 
version is believed to confine and break its discrete anomaly-free chiral symmetryPH 
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One can most easily see that this must happen by recalling the axial anomaly. 
As is well known, in quantum field theory (QFT) with Dirac spinors, even if the 
action is invariant under the axial rotation: 

V^^e iQ75 V^e 1 " 75 ^, (8) 

the partition function is, in general, not invariant. One way to understand this is due 
to Fujikawa^, who pointed out that although naively such a field rotation should 
certainly leave the fermion measure unchanged (as one expects that tr a(x)75 = 
for arbitrary a(x)), because the fermion phase space is infinitely dimensional, d = oo 
in our notation, it needs to be regularized and there exists no method of doing so 
respecting both the vector and the axial symmetry simultaneously. 

Trying to repeat the same argument on the lattice immediately leads to prob- 
lems. On a finite lattice, the theory is perfectly regularized and d is a finite integer. 
We apparently find a regularization scheme of the measure, which is invariant under 
both vector and axial rotations, leaving absolutely no room for an anomaly to arise. 
Applying this logic to gauge symmetries and gauge anomalies in chiral theories, we 
would similarly conclude that gauge anomalies do not exist on the lattice. Since in 
the continuum left- and right-handed Weyl fermions contribute to the anomaly with 
opposite signs and cancel if their charges are equal, the absence of anomalies on the 
lattice is only possible if every Weyl fermion that survives in the continuum limit is 
accompanied by a doubler with precisely the same charge but of opposite chirality. 
In other words, the fermion species are automatically doubled. The impossibility of 
preserving an exact chiral symmetry of the form ^ on the lattice without fermion 
doubling is the (in)famous fermion doubling problem. 

2.4. Ginsparg- Wilson fermions: an elegant solution 

According to the above discussion, it should be clear that the only way to avoid the 
fermion doubling problem on the lattice is to explicitly break the axial (or chiral) 
symmetry slightly so that the action is only approximately invariant. One hopes 
that in the continuum limit, this approximate symmetry can be fine-tuned to be 
exact and that the anomaly is properly reproduced as a remnant of the explicit 
breaking at the cutoff scale, see^H Methods of this kind exist and are widely used 
in vector-like theories, e.g. QCD. 

However, having only an approximate symmetry in lattice QFT has many dis- 
advantages. For our purposes, the most important is that it forbids pure chiral 
actions — or any action with a vectorlike fermion content (meaning equal number 
of left- and right-handed Weyl fermions with the same charge so they form vector 
multiplets), but consisting of two separated chiral sectors. Because the chiral sym- 
metry is explicitly broken by the regulator, fermions with opposite chiralities must 
necessarily couple. 

Ginsparg and Wilson (GW) proposed^ an elegant alternative in 1982. They 
suggested replacing the Dirac operator on the lattice by the so called GW operator, 
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which satisfies the following two conditions: 

{D, 75 } = aD l5 D, ( 75 £>)t = l5 D. (9) 

Here, a is the lattice spacing. In the continuum limit a — > 0, the rhs of the above 
anti-commutator vanishes and D approaches the usual Dirac operator. 

Using the Ginsparg- Wilson operator in the fermion kinetic terms is sufficient 
to eliminate the fermion doubling problem. It works, of course, because D explic- 
itly breaks the chiral symmetry (as it does not anticommute with 75). What is 
really beautiful, however, is that the GW operator allows to define an exact chiral 
symmetry on the lattice. To see this, define a new 75 matrix by: 

75 = (1 - aD) 75 , (10) 

and observe that (|9]) implies that the following two equations are exact: 

75^ + ^75 = 0, f 5 =l. (11) 

The first equation says that actions like S — ipDip are invariant under the GW axial 
rotation: 

^V^ 5 , V^e*" 75 ^, (12) 

which is an exact symmetry of the action at a finite lattice spacing. In the continuum 
limit when the lattice spacing vanishes, 75 — > 75 and the traditional axial rotation 
Q is recovered. 

We note that while the GW relation was proposed in 1982, no explicit form of 
an operator obeying ^ was known until 1997, when it was found by Neuberger 
as a result of a remarkable development) 41 ! 42 ! 43 ! 44 ! 45 ! We will not need the explicit 
form of the GW operator in this review, but only the properties (9p2 ) it satisfies; 
see^ for a recent introduction. 

Now, 7! = 1 implies that v B define two projection operators, and therefore 
one can build, out of Dirac spinors, the "GW chiral fermions" as follows: 

Vl/r = V 2 , Wl/r = 2 V- (13) 

This definition for the spinor tp^ /r is identical to the usual one in the continuum, but 
for ip, it is slightly different — but the difference disappears in the continuum limit. 
Thus, one can define "chiral" theories on a finite lattice, free of fermion doubling, 
using these GW chiral fermions. For example, a single massless free fermion can be 



described by the action S — ip L DipL, where ip L , defined by (13), is not the usual 
Weyl fermion that we are familiar with in the continuum. If there are no gauge fields, 
the measure can also be easily defined, as we will discuss later. The difficulties that 
arise when gauge fields are present are discussed in Section |2.5| below. 

Let us now explain how the correct chiral anomaly is derived in the GW for- 
malism. In a sense, the anomaly works out just as Fujikawa pointed out in the 
continuum: the anomaly exists because the fermion measure is not invariant. To see 
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how the index theorem works out in vectorlike theories, let us again start with the 
partition function for a single Dirac fermion: 

ThpVi/j exp (ipDip) , (14) 

where ip and ip are both Dirac spinors. In any vectorlike theory, the fermion mea- 
sure in the path integral is trivially defined on lattice. One simply chooses any or- 
thonormal basis of convenience (momentum eigenstates or lattice-site eigenstates, 
for example), expands the spinors with the chosen basis vectors with Grassmann 
expansion coefficients, and defines the fermion integral by the usual Grassmann in- 
tegral. When the lattice has a finite size, such an orthonormal basis is always finite 



and even dimensional, with dimension we denote as d. The operator D in ( 14 1 is 
the GW operator introduced before, satisfying the conditions From now on, 
we implicitly understand that some background gauge field is present and the GW 
operator D is a functional of the gauge field configuration, because it must approach 
the covariant Dirac operator in the continuum limit. As already mentioned, the ex- 
plicit form of D can be complicated and is not needed for the general discussion 
here. We also set the lattice spacing a = 1 from now on. 

In the four-dimensional Weyl representation of 7-matrices, 75 is given by: 

(75)2^ = ( 1 1 ) , £ is any lattice site, (15) 



and ('J5)xy — if x 7^ y. On the other hand, 75 is defined through equation ( |10[ ). 
Again, it suffices for now to recall that 75^ = 75 and -fs 2 = 1 (following from the 
GW relations) and that 75 — > 75 in the continuum limit. 



Classically, the theory ( 14 1 has two global symmetries. Under the vector rotation, 
the spinors transform as: 

^e iQ V, i>^i>e~ la , (16) 
and under the axial rotation, they transform as: 

4> -> e lai5 ip, i> -> ^e iai ~ 5 , (17) 

where a is some space-time constant. It is well-known that the axial symmetry 
suffers from a quantum anomaly when the path integral is properly regularized in 
the continuum. On the lattice, this anomaly appears precisely due to the fact that 
the axial rotation performed above does not leave the fermion measure invariant. 
Suppose (u±, 1*2, • • • , u<i) forms an orthonormal spinor basis of dimension d, and 
expand the fermion fields as follows: 

d d 

j> = y^ J a n u n , ip = b n ul . (18) 

n— 1 n—1 

The fermion measure is then just given by: 

d d 

V^Vi> = Yl db n Yl da n . (19) 
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Under the axial rotation of eqn. ( [17] ): 

ip' = ^2 a 'n u n = e M75 ^2 a ™ w ™> ^ = = X! ^« u » e 



(20) 



we find that to leading order in a: 

a' n = a n + iay~* y a m (^, 75U m ), b n = b n + ia'^b m {u\ n i b ,u n ). (21) 



Therefore, the Jacobian of the infinitesimal field transformation (20) is given by 



J- 



1 + ia^2(uli 5 ,u n ) 



la 



(22) 

=(1 + iaTrj^) • (1 + iaTi'75) = 1 + iaTrfs. 

We have used the fact that Tr75 = 0. The trace of 75, on the other hand, does not 
vanish in topologically nontrivial gauge-field backgrounds. The properties 75^ = 
75 and 75 2 = 1 alone are sufficient to deduce that the eigenvalues of 75 are ±1. 



Therefore, the trace appearing in (22), 



Trf 5 = n+ - n_, (23) 

must be an integeiQ Here n± are the dimensions of the invariant vector spaces 
spanned by the eigenvectors of 75 of ±1 eigenvalues respectively. Obviously n + + 
rt_ = d. This is the lattice version of the index theorem in the GW formalism. 

There is a significant amount of literaturd 46 ! 47 ! 48 ! 49 ! 50 ! on many aspects of the 
expression ( |23[ ) that will not be discussed here. For now, we only mention that one 
finds, using the definition of 75 and the explicit expression for D, expanded in the 
continuum limit 

T{K) 



tr(75xa 



32tt 



e »vapF£ v {x) + higher dimensional terms . 



(24) 



for 4d Dirac fermions transforming in a representation 1Z of SU(N), with quadratic 
Casimir T(1Z) (1/2 for the fundamental), or: 

1 



tr(7 5:Ea .) = — — envF^ u (x) + higher dimensional terms . . . , (25) 

27T 

for a 2d Dirac fermion of unit charge under a U(l) gauge group; in both expressions, 



x denotes an arbitrary lattice site. The expressions (24), (25) show that the usual 
index theorem of the continuum is recovered. 

To summarize this Section, the advantage of the GW formalism is that it elim- 
inates the fermion doubling problem and at the same time defines an exact chiral 
symmetry on the lattice, which reduces to the chiral symmetry ^ in the continuum 
limit. Consequently, fermions with opposite GW chiralities can be disentangled in 
the action, allowing the definition of "chiral partition functions," which play an 
important role in our further discussions. 



'It must be an even integer since n_|- + n_ = d is even. 

J We use "Tr" to denote trace over spacetime, gauge, and representation indices, while "tr" only 
includes spinor and representation trace. 
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2.5. Problems with Ginsparg- Wilson fermions in chiral gauge 
theories 

Sadly, the beauty of the GW formalism is not the end of the story of defining chiral 
gauge theories and troubles arise the moment we consider dynamical gauge fields. 
The discussion of Section [2.1| about the phase ambiguity of chiral partition functions 
finally becomes relevant. Again, consider the partition function: 

Z = J D^} L DipL exp (V^l) , (26) 



but let now ip L be the GW chiral fermions, as defined in (13). To define Z, one 



chooses sets of orthonormal basis for both ipL and ij) L , denoted as \uj \ and {vi}, 
respectively, and then evaluates Z as det(i>]D«i) (see also Section 



4.1.2 



for more 

details). As discussed in Section 2.1 the freedom of rotating either of the two basis 
by a unitary matrix renders the phase of Z ambiguous. 

What is changed in the GW formalism is that with dynamical gauge fields, one 
no longer has the option to just fix a set of vi and then stick to them. This occurs 
because the definition of i\) L involves 75, which depends on the Ginsparg- Wilson 
operator D, which in turn depends on the gauge field background. When the gauge 
field varies, the subspace in which ip L lives rotates along. The same set of v^s 
inevitably fails to span it. It becomes mandatory to choose a different orthonormal 
basis {vi} for every gauge field configuration, and all of them are subject to arbitrary 
unitary rotations. One finds, in the Ginsparg- Wilson case, not a single ambiguous 
phase but a t/(l) bundle fibered over the entire gauge field configuration space, or 
more rigorously speaking, a U(l) bundle fibered over the gauge field configuration 
space modulo all the gauge transformation — since the phase of Z should be at least 
gauge invariant (for an anomaly-free theory) to be considered interesting. Finding a 
particular way of fixing the phase of Z for all gauge field backgrounds is equivalent to 
finding a global section of this U(l) bundle. It is not obvious whether such a global 
section exists and whether it is unique when it does. The problem of determining 
the phase of Z is often referred to as the "fermion measure problem" in lattice chiral 
gauge theories. 

Although there is no natural choice of a connection on this £7(1) bundle, there 
exists a natural (i.e., independent of the basis vectors) definition of its curvature, 
solely determined in terms of the GW operator D. This shows that the phase of the 
partition function can not be chosen at random — in particular, if the phase is taken 



to be gauge-field independent, the partition function (261 is not even gauge invari- 
ant, even in the anomaly- free case|^] Integrating this curvature over any nontrivial 
cycles in the gauge field configuration space, Neuberger^ and LiischeJ ^ 2 ! 53 ! 54 ! found 



k Disc ussing this t opic in detail is beyond the scope of this review. We refer the reader to the original 
paperi 51 52 53 54 or the reviewiHIHIHl For an illustration (within the Wilson-line subspace of the 
gauge-field configuration space) of the topological obstruction to defining Z in the anomalous case 
and an explicit construction of the phase of Z in the anomaly-free case, see^l 
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that it vanishes precisely when the gauge-anomaly cancellation condition is obeyed. 
This result indicates that in the anomaly-free case a global section can be found. A 
proof that a global section — one that ensures gauge invariance, locality, and smooth- 
ness of Z as a function of the gauge background — can be found and is unique, up 
to usual irrelevant counterterm ambiguities, in the case of an anomaly-free abelian 
gauge group, was given irP^. Despite Liischer's existence proof, a practical way of 
implementing the theory on a computer is still missing; however, see the work of 
refsl 55 ! 56 ! on U{1) and SU(2) x U(l) chiral gauge theories. At the moment of writ- 
ing, a generalization of the proof to general anomaly-free nonabelian chiral gauge 
theories is not known. 

Apart from the fermion measure problem, there are other issues with the formu- 
lation of lattice chiral gauge theories with Ginsparg- Wilson fermions by the path 



integral (26): 



(1) The issue of summation over different topological sectors. This issue arises be- 
cause in different topological backgrounds the dimensionality of the fermion 



path integral in (26 1 changes. This is because the dimensionality of the ip L 
space is equal to trP + , which has different values in different topological sectors 
(on the other hand, the dimensionality of tpL space is gauge-background inde- 
pendent, equal to trP+). What is not known is how to determine the relative 
normalization of the contribution of different topological sectors to expecta- 



tion values of observables via (26 1. This is not an issue in vectorlike theories, 
where the dimension of the iPl + r space is trP|_+trP_ = trP++trP_ and is thus 
equal to the dimensionality of the ipL+R space for any gauge background. In 
this context, we note the proposal oP^ to determine the absolute value of the 
normalization factors, but not their phase, starting from a real-representation 
vectorlike theory (however, it is not clear if this normalization is consistent 
with cluster decomposition). If the "decoupling of mirror fermions" approach 
discussed in this review is successful, this normalization issue is resolved auto- 
matically, since the measure is that of the vectorlike theory. 
(2) The unwanted- CP- violation problem. Under the condition of locality, the gauge- 
field dependent projectors P± have an asymmetry between fermions and anti- 
fermions. As this asymmetry provides the origin of fermion number violation 
in the chiral theory, it is a welcome feature. However, at the same time it 
breaks CFl 58 l 59 l 6 H This CP violation afflicts high-momentum modes and might 
be (naively) expected to disappear in the continuum limit; furthermore, it may 
also show up in instanton-induced fermion-number violating amplitudes^. We 
stress that the effect of this CP violation in the quantum continuum limit is 
not completely understood and that the CP violation issue also afflicts the 
mirror-decoupling approach (if standard projectors and CP are used), where 
its implications await further study. Note that the above remarks refer to the 
continuum CP transformation: recently refP^ argued that in the trivial topolog- 
ical sector it is possible to define a lattice-modified version of CP, (which, like 
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the modified chiral symmetry, reduces to the usual continuum CP-transform), 
under which the chiral action (26) (and measure) is invariant; see alsd 62 l 63 l 



3. Chiral from vectorlike? 

3.1. The idea of decoupling the mirror fermions 

Our approach to tackle the long-standing problem of defining lattice chiral gauge 
theories is to try to work around the toughest part. Realizing that defining the mea- 
sure of an arbitrary chiral gauge theory explicitly is beyond our reach, we propose 
to start with a vectorlike theory insteadj^The measure and the action of a vectorlike 
theory can be unambiguously defined. 

Thus, the vectorlike theory we consider should be described by an action S which 
splits into two separate chiral parts: 

S = Slight + S m i rror . (27) 

The fields that participate in the two parts of the action will be called the "light" 
and "mirror" fields, respectively. The light sector describes a chiral gauge theory 
that is the target theory we wish to obtain in the end, while the mirror sector 
contains all the fermions with the wrong chiralities. As explained earlier, such a 
clear separation is only possible on the lattice using the GW formalism; hence, the 
so-called chiral fermions refer to the GW chiral fermions. 

The hope is to design the mirror sector in such a way that all the mirror fermions 
of wrong chiralities contained there become heavy and disappear in the IR limit. If 
the gauge symmetry remains unbroken by the mirror dynamics, we would obtain 
a chiral gauge theory described by the action Slight- Since the measure is that of 
the vectorlike theory, we never have to worry about choosing the ambiguous phase, 
thus circumventing the difficulty of defining the chiral fermion measure explicitly. 
We call this idea "decoupling of the mirror fermions." 

To better explain what we are proposing, let us first work in continuum terms 
and explain the desired features of the mirror dynamics. We will use the chiral 
SU(5) theory already discussed in the Introduction to illustrate the idea. We use 
two-component Weyl-spinor notation to describe the desired "light" fermions of the 
target SU(5) theory: 

light : ifc ~ 5*, x O a ~10, Ca ~ 1 , (28) 
where i denotes an S/7(5) (anti-)fundamental index and a = 1,2 — an SL{2, C) 



index. We also introduce the "mirror" partners of ( 28 1 



mirror : r) ia ~ 5, p» ~ 10* , ~ 1* . (29) 

In this notation a Dirac mass term for the 5 + 5* would be of the form rmp la r]i a + 
h.c. We also introduced a gauge-singlet Dirac fermion, with Weyl components 

'We should note right away that this is not a completely new idea. However, when coupled with 
exact lattice chirality, it acquires many new aspects, see Section 13.41 for discussion and references. 
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This is a field whose £ component will play an important role in the strong mirror 
dynamics (its SU{b) representation was denoted by 1* simply to distinguish from its 
partner) [^ As already stated, we imagine that on the lattice the chiral components 



are defined as appropriate for GW fcrmions ( 13 ) 



Now, to decouple the mirrors (291, we add non-gauge interactions involving only 
the mirror fields. Since the mirror fermions are in a chiral representation, there are 
no mirror-fermion gauge-invariant bilinears that one can write down. However, three 
gauge-invariant four-fermion interactions can be written down. In group theory 
terms, these invariants are 5-10*-10*-10*, 1-5-5-10*, and 1-1-1-1. More explicitly, 
the mirror action is: 

Smirror = ^ff^pP ^P^Jklmn + WVa irf P% + ^"SettfZfi + h.C, (30) 

where we omitted the mirror kinetic terms, which have the standard form. The rea- 
son for introducing the mirror interactions is to decouple the mirror fermions 
while preserving the gauged SU (5) symmetry. This might be possible if the four- 



fermi interactions (30), when taken strong, lead to the formation of SU(5) invariant 
(or vectorlike) mirror composite states, which can acquire mass without breaking 
SU(5). For example, a possible composite of the mirror fermions is the rjijp (5-5- 



10*) invariant appearing in (30). This composite state can acquire a large Dirac 
mass by pairing with the singlet mirror field £ and can thus decouple from the 
low-energy physics without breaking the SU(5) symmetry. Similarly, at strong cou- 
pling, all mirror fermions are expected to be bound in massive singlet or vectorlikc 
composites and decouple from the infrared dynamics. 

Thus, the desired spectrum of the theory consists of heavy mirror states, with 



mass of order of the ultraviolet cutoff (see Section 3.2 1. The SU (5) gauge interactions 
are only a spectator to the strong mirror dynamics. The gauge interactions couple 
to both the light and mirror states, but the gauge coupling at the cutoff scale is 
expected to be weak and is not expected to cause the heavy mirror states to become 
light or the massless light states to become heavy (the latter are protected by their 
unbroken exact chiral symmetries). Thus, at scales below the cutoff, the spectrum 
and the global chiral symmetries of the light fermions would be exactly the ones of 
the desired chiral gauge theory. 

There are several pieces of evidence that the previous two paragraphs are not 
complete fantasy. 

The first — and most pertinent to our approach — is the existence of strong- 
coupling symmetric phases of lattice four-fermi or Yukawa models. In the following 
Section [3. 2| we consider a simple toy example of a strong-coupling symmetric phase 
in a theory with naive lattice fermions. In this phase, the spectrum has only mas- 
sive states and an unbroken global symmetry. The nature of similar strong-coupling 
symmetric phases in theories with GW fermions is not yet completely understood 

m An entire gauge-singlet Dirac multiplet was added to make sure the fermion representation ) |28| 
[291 is vectorlike (even though this is not strictly necessary in the singlet sector). 
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and is the subject of later Sections. 

The second — less directly relevant to our approach, but nonetheless tantalizing — 
motivation is provided by the Af = 1 supersymmetric examples where chiral and vec- 
torlike gauge theories are related by Seiberg dualities. There, a theory with charged 
matter content in vectorlike representations — but with some chiral couplings to sin- 
glets in the superpotential, not unlike (30 1 — is argued to have an equivalent infrared 
description in terms of a theory with chiral matter content. This example shows 
that, at least in the supersymmetric "theory space," gauge theories with chiral 
and vectorlike matter content may be related. The first example of chiral-nonchiral 
Seiberg duality was found by PouliolpS and was later generalized by Pouliot and 
Strassler^S. They considered the vectorlike J\f = 1 supersymmetric SO(8) gauge 
theory with one chiral superfield in the spinor and N chiral superfields in the vector 
representation (both are real, eight-dimensional representations), along with some 
gauge-singlet chiral superfields with chiral Yukawa couplings. They argued that this 
vectorlike theory has a dual infrared description in terms of a chiral SU (N) theory 
(6 < N < 16) with a symmetric tensor and N + 4 antifundamcntal chiral multiplets 
without any superpotential. The two dual theories — the chiral and the vectorlikc 
one — flow to an interacting conformal field theory. The chiral primary operators of 
the two theories and their scaling dimensions (which can be calculated in super- 
symmetry even at strong coupling) can be mapped to each other and the duality 
map can be shown to pass some rather intricate consistency checks. Furthermore, 
as in Seiberg duality between vectorlike theories, a free magnetic phas^] occurs for 
N > 17 — where the dual vectorlike SO(8) theory loses asymptotic freedom. Thus, 
the asymptotically free chiral SU(N — 4) theory flows in the IR to a theory of free 
50(8) vectorlikc quarks and gluons (and some gauge singlets), which can be viewed 
as the composite quasi-particles of the strongly-coupled chiral theoryj^] 



3.2. A toy model of decoupling in strong- coupling symmetric 
phases. 

In this Section, we consider a simple toy model illustrating the idea of lattice strong- 
coupling phases in theories with multi-fermion interactions. The simple model be- 
low was invented solely for the purpose to give a short introduction to strong- 
coupling symmetric phases in multi-fermion or Yukawa theories on the lattice. 
Se€ )Bb|S7|B8|S9|7U|7i|72|73|74| for detai i ed stu dies of related models, where the entire 



n See the already mentioned^ for references. 

° It would be interesting to know if examples of chiral-nonchiral dual pairs with a free magnetic 
phase described by a chiral theory exist (this question is not easy to answer, as there is no known 
"algorithm" for finding Seiberg duals). If so, the chiral theory would arise as an effective weakly- 
coupled IR description of the vectorlike theory: a more ambitious goal than what we are trying to 
accomplish with the "mirror decoupling" idea (as our aim is to only make the mirror fermions heavy 
without breaking the gauge symmetry, which is only a spectator to the strong mirror dynamics). 
In a chiral magnetic dual, the gauge fields and the fermions of the chiral theory would appear as 
composite objects of the vectorlike theory — and would thus offer an sintriguing possibility. 
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phase diagram is also discussed; for brevity, we only focus on the strong-coupling 
symmetric phase as it is the one relevant for us. 

We will consider our lattice toy model in the Hamiltonian formulation (where 
space is a lattice, but time is continuous). The fields are one-component fermion 
fields tl> a {x), a — 1, ... 4, living on the sites x of a lattice of any dimensionality. They 
obey the canonical anticommutation relations: 

{ip a (x),^l(x)} = S ab S xy , (31) 

and have a local 4-fermion interaction: 

Hint = £ A {Mx)Mx)Mx)Mxy bcd + h-c) • (32) 

There are also some hopping terms in the Hamiltonian, whose particular form is 
not important in what follows, say: 

H km = ^i>l(x + MM*) +••■, (33) 

where /t is a unit lattice vector in the fi-th direction. The lattice spacing has been 
set to unity. The Hamiltonian preserves an SU(4) global symmetry The analogy 



we want to keep in mind is that the 4-fermion Hamiltonian ( 32 ) is similar to the 



mirror interactions in the SU(5) theory (see eqn. (30)) and that the global SU(4) 
symmetry is the analog of the to-be-gauged SU(5). 

We will study the limit A > 1, when the four-fermi coupling is strong in lattice 
units. In this limit, the Hamiltonian of our toy model can be easily diagonalized, 



since to first approximation one can neglect the hopping term ( 33 ) , decoupling all 



lattice sites in (32). Thus, it is possible to find the eigenstates of the Hamiltonian 
( 32 1 on every lattice cite and then build the ground state of the theory as a direct 
product of all single-site ground states. To diagonalize the single-site Hamiltonian, 
we denote ipb — ► a>b, i>\ - > a \, where oj,, a\ are four fermion creation and annihilation 
operators obeying the usual canonical anticommutation relations {afc,aj} = #& c . 
The single-site Hamiltonian is a simple quantum mechanics problem of four fermions 
with an interaction that preserves fermion number F modulo 4 and an SU(A) global 
symmetry: 

H a = A (a a a b a c a d + a\a\a\aty e abcd . (34) 

The Hilbert space of the single-site problem has 2 4 states which can be decomposed 
as SU(A) representations of given fermion number F, with |0) denoting the Fock 
vacuum, as shown in Table [T] Conservation of F(mod)4 implies that only H only 
has nonzero matrix elements between the states |1) and Thus, we find that the 
state |1) — |1') has energy —A, the state |1) + |1') has energy A, while the states 
|4), |6), and |4*) have zero energy. Hence, in the A — > oo limit, the ground state of 
the single-site Hamiltonian is unique and is an SU(A) singlet. The energy levels of 
H are shown on Fig. 1. 
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Table 1. The 2 4 states of the 
single-site four-fermi Hamilto- 
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Fig. 1. Spectrum of the single-site hamiltonian ( |')4| of the A — > oo limit our toy model | |32| |33[ . 
The ground state is unique, has a gap A in lattice units, and preserves 5(7(4). Taking hopping 
terms into account causes massive excitations to propagate between adjacent sites on the lattice. 

We can now build the A — > oo ground state of the entire lattice theory by simply 
taking the product of the |1) — |1') ground states of every lattice site. This state 
is unique and is an SU (4)-singlet. The name "strong-coupling symmetric phase" 
should be clear now: the ground state of the four-fermi theory at A — > oo is an 
SU(4:) singlet and all excitations are gapped, with mass gap proportional to A in 
lattice units. Taking hopping into account can be organized in a 1/A strong-coupling 
expansion; it causes the localized states to propagate from site to site, leading to 
propagating states of mass ~ A. Since the strong-coupling expansion has a finite 
radius of convergence, the picture of an SU (4)-singlet ground state and only massive 
excitations is a true representation of the ground state of the finite- A theory, for 
sufficiently large A> 1. A transition to a broken phase is expected to occur at some 
critical value A c ~ 1; this has been studied in the literature quoted in the beginning 
of this Section. 

The strong-coupling symmetric phase in this toy model has features similar to 



what we desire of the mirror dynamics ( 30 ) : for sufficiently strong mirror coupling it 



gives rise to only heavy mirror states and the to-be-gauged global SU (4) symmetry 
is unbroken. The physics of the strong-coupling symmetric phase of this toy model 
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is thus "trivial," as there are no states lighter than the ultraviolet cutoff and thus 
no nontrivial infrared physics. Upon gauging the SU{A) symmetry, however, there 
would be nontrivial infrared dynamics — that of the unbroken pure gauge theorjjp] 

In his context, we note that the question of the "triviality" of four-fermi (or 
Yukawa) interactions in 4d is often raised when their strong-coupling symmetric 
phases are mentioned. In our context, triviality — interpreted as the absence of any 
long-distance physics in the strong-coupling symmetric phase — is to be advanta- 
geously exploited, since the goal of the strong mirror interactions interactions is to 
keep all mirror fermions at the cutoff scale. 

For the mirror-decoupling approach to chiral lattice gauge theories, the most 
important and still not completely answered question is whether models with lat- 
tice four-fermi mirror interactions formulated with Ginsparg- Wilson fermions have 
strong-coupling symmetric phases which share the same desired properties: an un- 
broken global symmetry and no massless mirror states. We already noted that using 
GW fermions allows us to separate the light and mirror modes in an unambiguous 
way while preserving all the chiral symmetries of the desired target theory. However, 
the dynamics of GW fermions is more complicated than that of the simple 4-fermi 



model with local interactions (32). In particular, a strong-coupling expansion that 
preserves the chiral symmetries appears to be difficult to perform, because of the 
exponential-onljCHZZl locality (which holds in the "admissible" part of the gauge- 
field configuration space) of the exactly chiral lattice fermions. Thus, the numerical 
study of the strong-coupling phases appears the only way to proceed. We will de- 
scribe the status of these studies in the following sections. 



3.3. Mirror global symmetries and 't Hooft anomaly matching. 

In this Section, we explain why a necessary condition for the decoupling of the 
mirror fermions is that the mirror interactions explicitly break all mirror global 
symmetries, including the ones that would have a gauge anomaly after gauging 
SU(5). To this end, note that in the limit when the SU(b) gauge coupling and the 



interactions in (30) are turned off, the mirror theory has an SU(5) x U(l) 3 global 
chiral symmetry. The U(l) 3 are phase transformations acting on each of the three 



mirror multiplets n,p,^. The interaction terms in (30) explicitly break [/(l) 3 , but 
preserve the to-be-gauged SU(5). 

The mirror theory is a theory with a chiral matter content and its unbroken 
chiral global symmetries have 't Hooft anomalies. The 't Hooft anomaly-matching 
argument is usually applied to asymptotically-free gauge interactions that become 
strong in the infrared. However, 't Hooft anomaly matching also holds when strong 



non-gauge interactions on the lattice are considered (such as the ones in (30)), 



when the theory is formulated with exactly chirally-symmetric GW fermions. As we 



p For a recent study of the decoupling of scalars in a strong-coupling symmetric phase in a simple 
toy model with gauge fields, 
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explain in detail in Sections 4.1.4 and |4.2| the reason is essentially that in theories 
of chiral fermion formulated via GW fermions, the anomaly is independent of the 
coupling strengths and the details of the lagrangian (see Section 4.1.2). 

For now, taking the anomaly-matching condition for granted, note that the 
SU(5) symmetry is respected by the strong mirror interactions in Smin-or, but since 
the mirror spectrum is anomaly-free, there is no SU(5) 3 't Hooft anomaly to satu- 
rate. However, if one of the couplings Ai.2,3 in ( 30 1 was set to zero, there would be an 
additional anomaly-free chiral U(l) respected by the mirror interactions. This extra 
U(l) would have some 't Hooft anomalies (U(l), U(l) 3 , or U(1)SU(5) 2 , depending 
on which coupling was set to zero). At zero mirror couplings, these anomalies are 
due to the elementary massless mirror fermions. The 't Hooft anomaly-matching 
conditions imply then that this anomaly has to be matched by the spectrum of the 
strongly-coupled mirror theory — and thus light mirror states would be guaranteed 
to exist (whether they are bosons or fermions depends on the dynamics and can 
not be inferred from anomaly-matching arguments alone). Thus, we conclude that 
symmetry reasons alone require that all mirror global symmetries — except the one 
to be gauged — have to be explicitly broken, to prevent the appearance of massless 
mirror states required by 't Hooft anomaly matching. 

We stress that the absence of exact global chiral symmetries in the mirror the- 
ory is only a necessary condition — it does not guarantee that the mirror spectrum 
is necessarily free of massless fermions. The absence of anomalies to match simply 
means that any massless mirror states would have to come in anomaly-free rep- 
resentations under the unbroken SU(5) — as would happen, for example, at weak 
mirror couplings Ai,2,3- However, our study of a two-dimensional model in Section 
|4.3| shows that, at strong mirror coupling, the mirror spectrum saturates anomaly 
matching with the smallest required number of massless fermions and that mass- 
less fermions in anomaly-free representations do not occur when all mirror global 
symmetries are broken. Needless to say, we take this as an encouraging sign — but 
by no means a proof — of the viability of this approach. 

To end this Section, we note that another way to see that mirror global sym- 
metries that would be anomalous after gauging 577(5) have to be explicitly broken, 
is to consider what would happen in a topologically nontrivial background, as dis- 
cussed already irP. In an instanton field, all chiral fermions, mirror and light alike, 
have zero modes as dictated by the index theorem. In the SU(5) theory that we use 
as an illustration, this is shown on Fig. 2. In order to lift the mirror-fermion zero 
modes, one must include a mirror interaction that explicitly breaks the anomalous 
mirror global chiral symmetry. In this case, the Ai term in ( 30 ) allows all mirror 
zero modes to be lifted and the target theory 't Hooft vertex to be reproduced. 

Note the analogy with the Standard Model: in QCD, quark masses explicitly 
break the anomalous chiral symmetries violated in SU (3)-instanton backgrounds 
and allow heavy quarks to completely decouple from the instanton vertex. In con- 
trast, in the electroweak theory, Yukawa couplings respect the anomalous B + L, 
violated in an SU{2)l instanton background, and thus heavy quarks do not decouple 
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Fig. 2. In order to lift the mirror-fermion zero modes in an instanton background, the mirror 
interactions have to break the anomalous mirror global symmetries as well. 

from the baryon-number violating 't Hooft vertex^. Thus, in order to be success- 
ful, our mirror-decoupling approach should resemble the case of QCD, where the 
anomalous chiral symmetry of the heavy quarks is explicitly broken. 

3.4. Comparison to "pre- Ginsparg- Wilson" attempts of 
mirror / doubler decoupling 

This analysis of the toy model of Section |3.2| can be also performed for realistic 
models in four spacetime dimensions (and can also be done in an Euclidean lattice 
formulation). The idea of the analysis is the same, but the calculations become more 
cumbersome in view of the proliferation of indices to keep track of. 

In fact, such a study was done by Eichten and PreskilP^ for the St/ (5) model 



of Section 3.1 They studied the St/(5) model with local four-fermion interaction 



terms of the form Ai5-10*-10*-10* +A 2 l-5-5-10* + h.c. (essentially the ones of 



( 30 1 with A3 = 0) . All fields were taken to be naive lattice Weyl fermions, in order 
to preserve the chiral symmetry. In this 1986 formulation it is the 15 doublers (on 
an Euclidean lattice) of the naive Weyl fermions that play the role of our "mirrors." 
Eichten and Preskill's ultimate goal — similar to ours — was to find a phase where 
the doublers decoupled and the light components remained exactly massless. 

The two couplings, Ai and A2, were motivated by: a.) having to lift the mirror- 
fermion instanton zero modes and b.) ensuring that the Euclidean path integral has 
a "static limit." The first condition was already explained in the previous Section 
3.3| The second condition was the reason to add the singlet fermion with coupling 
A2 — if it is absent, the Euclidean path integral over fermions with kinetic terms 
put to zero vanishes, indicating that no static limit and no sensible strong-coupling 
expansion exist (in a Hamiltonian formulation this would mean that the ground 
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state at each site is not unique — indicating that there would be no stability with 
respect to hopping perturbations). 

By an analysis which is ideologically identical to the one of Section 3.2 Eichtcn 
and Preskill showed that when Ai and A2 were taken strong in lattice units, 
the ST/ (5) global symmetry remained unbroken, but all components of the Weyl 
fermions became heavy. Thus, while a strong-coupling symmetric phase exists, both 
the light fermions and the doublers are heavy, absent any additional couplings. They 
proceeded by adding two more four- Fermi terms with couplings n , T2- These terms 
have the same form as in ( 30 ) but now include lattice derivatives designed to lift the 
light-doubler degeneracy. The hope was to find a region in the enlarged 
parameter space of the model, where there would be a strong-coupling symmetric 
phase — but now with doublers having a cutoff-scale mass and exactly massless light 
fields in a 5 + 10* representation of the unbroken SU(5). The model with these 
extra couplings was analyzed irPHl, using all available analytic techniques. The con- 
clusion was that in all accessible corners of the phase diagram with unbroken SC/(5), 
either all fermions (light and doublers) were massless or all fermions obtained mass. 

While the study of refP^can not be considered as having given a proof, we believe 
that the apparent failure of this approach can be traced back to the absence of 
lattice chiral symmetries. To this end, note that in the Eichten-Preskill model, both 
the "light" and "mirror" (there: doubler) fermion components of the lattice Weyl 
fermions participate in the strong mirror interactions (recall that, in 1986, there 
was no known way to separate the chiral components from a lattice fermion field). 
While a mass splitting of the "light" and "mirror" (doubler) fermions occurs through 
the 7*1,2 derivative four-Fermi interactions mentioned above, exactly massless light 
fermions and heavy mirrors could at best be expected to occur for special, finely 
tuned, values of the "four-fermi Wilson" couplings r*j, as there was no lattice chiral 
symmetry protecting the light fields. A chiral symmetry could only be expected to 
emerge at certain values of rj . 

We stress that, in addition — as well as prior, see^for references — to the worlPSI 
there were many studies of mirror/ doubler fermion decoupling via strong Yukawa or 
four-Fermi interactions' 82 ! 83 ! 84 ! 85 ! 86 ! 87 ] (strong Yukawa and strong four-Fermi inter- 
actions can be mapped to each other; this was part of the analysis oP^l), reaching 
similar conclusions. The absence of lattice chiral symmetries was a feature of all 
these attempts. The lack of an exact chiral symmetry without doublers makes it 
impossible to split a lattice Dirac fermion field into chiral components (or separate 
a naive Weyl fermion field into light and doubler parts). As a result, any strong 
non-gauge interaction inevitably couples to both the "mirror/doubler" and "light" 
components of the fermion. The analysis of the strong dynamics — which is clearly 
not purely "mirror" anymore — of the chiral symmetries and their realization then 
becomes ambiguous. Most importantly, these formulations lacked a manifest sym- 
metry explaining why some chiral fermions should stay exactly massless while others 
could obtain mass; instead, chiral symmetries distinguishing "light" from "mirror" 
modes and protecting the "light" modes were expected to somehow emerge at spe- 
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cial values of the couplings (Aj. 2/^1,2 in the above example); see^ and references 
therein for related studies in Yukawa models. 

In contrast to the earlier approaches, the advent of exact lattice chiral symmetry 
in 1997 allows the construction of lattice actions where only the mirror fields partic- 
ipate in the strong interaction (a proposal along these lines was made^ in the finite 
domain- wall context, shortly before the advent of exact chiral lattice symmetry). 
The chiral symmetries of the mirror theory, their 't Hooft anomalies, and their real- 
ization can now be unambiguously defined and studied. The unbroken global chiral 
symmetries are exactly those of the target continuum chiral gauge theory and are 
manifest at finite lattice spacing and volume (as argued irP^) , leading to an elegant 
formulation of the mirror-decoupling idea. 



3.5. The big picture 

It is now straightforward to formulate mirror decoupling via Ginsparg- Wilson 



fcrmions, for example for the SU(5) model considered in Section 3.1 The action 



consists of kinetic terms for the vectorlike fermions transforming in the 5* and 10 



representation, the singlet fermion, plus the mirror action ( 30 1 . A light- mirror split 
of the lattice fermion action: 

S fermion — Slight ~t~ S m i rr0 r , (*^) 

is possible because of the existence of lattice chirality, which allows to split a vector- 
like fermion in chiral components at finite lattice spacing. The aim of the discussion 
here is to pose the various questions — answered or not — that arise in this approach 
for general models. 

The fermion partition function in a given gauge background can be written as 
follows: 

Z vector [A}= J dil>d$e? u ' M W +s ^""M , (36) 
where dtpdip is the well-defined fermion measure of the vectorlike theory, see 



eqn. (19 1. Exact lattice chirality allows (see Section 4.1.2) the partition function 



Zvector to be split, schematically: 

^vector 

[A] = Z Ught [A] x 

^mirror 

[A] (37) 

into light and mirror partition functions. The light partition function is that of 
the chiral 5* and 10 gauge theory, while the mirror partition function includes the 



mirror fermions plus their interactions. An explicit form of (37) in the 2d models 



whose dynamics we have investigated will be given Section |4.1.2| 



The split form of the partition function ( 37 1 is useful, because it allows us to focus 
attention on the non-gauge mirror dynamics described by Z mirror [A] , first without 
gauge fields (or with only perturbative gauge fields turned on). An investigation 
of the mirror dynamics without dynamical gauge fields is a sensible first step in 
understanding the possibility of mirror-fermion decoupling. 
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3.5.1. Questions about the mirror dynamics 

As already stressed, the question of defining the fermion measure, and thus of the 
chiral theory partition function, does not appear here. However, the price to pay is 
the introduction of strong mirror dynamics, about which many questions arise: 

(1) The multi-fermi (or Yukawa) interactions in the mirror theory formulated with 
GW fermions are only exponentially local, subject to the admissibility condition 
on the gauge field background 76 ! 77 ! Do strong-coupling symmetric phases exist? 

(2) In typical models, there are many multi-fermi/ Yukawa couplings in the mirror 



theory, e.g. ( 30 1 , needed to break all necessary mirror global symmetries. Is there 
a nontrivial phase structure as the ratios of these strong mirror couplings are 
varied? Is fine-tuning necessary to reach the strong-coupling symmetric phase? 

(3) Suppose one attempts to decouple an anomalous representation. Does the strong 
mirror dynamics avoid decoupling an anomalous representation in a symmet- 
ric phase (and thus creating an inconsistent anomalous long-distance theory)? 
How? 

(4) The Yukawa/multi-fcrmi couplings in S m i rror formulated with GW fermions 
define a complex partition function. Is the long-distance theory unitary? 

(5) Do mirror fermions in anomaly-free representations decouple in the symmetric 
phase at strong coupling, as suggested by symmetry arguments? 

(6) If the mirror dynamics works as hoped for, will this construction be useful for 
simulating any anomaly-free chiral gauge theory, in either 2d or 4 d? 



3.5.2. Answers — known and unknown 

We do not yet know the answers to all these questions. The answers below are based 
on symmetry arguments and evidence from studies of specific models: 

(1) Strong-coupling symmetric phases with GW fermions have been shown to exist 
in botrP2d ancEMUd. 

(2) There is a nontrivial phase structure in the strong-coupling regime^ in the 2d 
model we have studied. No fine-tuning is needed to reach the strong-coupling 
symmetric phase, as it exists for a wide range of couplings. 

(3) An anomalous representation does not decouple in the strong-coupling symmet- 
ric phase. A charged mirror fermion remains massless, realizing the minimal 
solution of the relevant 't Hooft anomaly matching conditional. 

(4) The result from (3) gives evidence for long-distance unitarity of the strong 
mirror interactions (as unitarity is crucial for 't Hooft anomaly matching). 

(5) We do not know the answer yet. However, numerical studies of mirror- 
decoupling of anomaly-free 2d models are feasible and should provide evidence 
for (or against) decoupling in the near future. 

(6) This will depend on the severity of the sign problem. Based on continuum 
intuition, a first test of the method with dynamical gauge fields is likely to be 
possible in 2d chiral U(l) or in the 4d j=3/2 Weyl fermion SU(2) theory. 
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In the rest of this article, we will explain in more detail our answers above. Most 
of them were obtained from numerical simulations. This is because of the lack of an 
obvious controlled expansion to study the GW-fermion strong mirror dynamics]^ 
Due to the high cost of simulations with exactly chiral fermions, such studies are in 
their infancy. On the other hand, since the issues considered are ones of principle, 
analyzing two-dimensional models is a sensible first step. Showing that the ideas 
work in 2d will not prove that 4d chiral gauge theories can be similarly formulated 
(this will require 4d mirror studies) , but the results are likely to provide insight into 
the relevant mirror dynamics. 

Another appropriate simplification is to neglect the gauge field fluctuations, 
since gauge fields play a spectator role to the strong mirror dynamics. If the strong 
non-gauge mirror dynamics gives the mirror fermions mass of order the lattice cut- 
off, in a manner similar to the one in the toy model, asymptotic freedom of the 
gauge interactions and the exact chiral symmetry of the light fermions lead us to 
expect that turning on dynamical gauge fields will not significantly affect the mirror 
dynamics or lift the masslcss modes. Questions about CP violation in topologically 
nontrivial backgrounds will have to be investigated, see Section |2.5| Including dy- 
namical gauge fields would only be warranted after the mirror decoupling in zero 
gauge background is demonstrated. 

4. Theoretical and Monte-Carlo studies of mirror decoupling via 
Ginsparg- Wilson fermions 

4.1. Strong- coupling phases and properties of chiral partition 
functions 

Most of the answers listed in Section 13.5.21 above were obtained in the framework 
of a toy 2d model, which we call the "1-0" model. This model is defined in Section 
14. UN below. 

We will also use the "1-0" model to illustrate various theoretical results, which 
are more generally valid: the splitting of a vectorlike partition function into chiral 
partition functions and an important result on the variation of chiral partition 
functions with the gauge background — the "splitting theorem" — a result without 
which many of the analytic and numerical studies would not be possible (see Section 
4X21. 



q Since, as already mentioned in Section [372] the GW-fermion multi-fermion/ Yukawa interactions 
are only exponentially local (subject to the admissibility conditiorEUZZl) . 

r The name "1-0" model is reminiscent of a general class of anomaly-free £/(l) chiral gauge theories 
in 2d— the "3-4-5", "1-1-1-2", etc., models. Recall that the anomaly cancellation condition for a 
U(l) chiral gauge theory in 2d states that the models should be "Pythagorean", i.e. the charges <j_ 
of the left-handed and g+ of the right-handed Weyl fermions should obey ^2 1— = Clearly 
a chiral gauge theory with "1-0" charges of the left and right moving fermions (the "light" or 
"mirror" sectors of the 1-0 model) would be anomalous. 
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4.1.1. The "1-0" model 

The "1-0" model is a Yukawa-Higgs-GW-fermion model: a U (1) two-dimensional lat- 
tice gauge theory with one Dirac fermion tp of charge 1 and a neutral spectator Dirac 
fcrmion x- Considering this theory was motivated by its simplicity: it is the minimal 
Higgs-Yukawa-GW-fermion model in two dimensions with chiral Yukawa couplings. 
It holds the promise to yield, at strong Yukawa coupling, a strong-coupling symmet- 
ric phase and is, at the same time, amenable to numerical simulations not requiring 
the use of extensive computing resources. 

The action of the "1-0" model, with the gauge kinetic term omitted, is: 

S Slight ~l~ S m i rror 

Sught = ~ (ip+ -Di'i>+)- (X- ■ A) • X-) 

Smirror = S K - (lp_ ■ D X ■ ijjj) - (x + ' D ■ X+) 

+ y{(^--<f>* ■x + ) + (x+-<t>-^-)} (38) 
+ yhi (VZ • 072 • x+) - (x+ ■ 72 • 4>* ■ } • 

The chirality component^] for the charged and neutral barred fermions are defined 
by projectors including the appropriate Neuberger-Dirac operators (charged Di and 
neutral D ), for example ip± = ^(1 T 7s)/2- The brackets denote summation over 
the lattice sites as well as a spinor inner product. The field cf> x — e lr):c , |?7 X | < 7r, is 
a unitary Higgs field of unit charge with the usual kinetic term: 

s - = \ E E t 2 - ( & u *>*+* ^+/» + h - c - )] . ( 39 ) 

x ft 

where U x ,x+)i = e lA ^ xS > . We will call the fermion fields that participate in the 
Yukawa interactions the "mirror" fields — these are the negative chirality component, 
ip_, of the charged ip, and the positive chirality component, x+> of the neutral x — 
while the fields ip+ and X- wm De termed "light." 



When the Yukawa couplings are absent, at y = 0, (38) is just the Schwinger 
model with a neutral Dirac spectator fermion and a unitary Higgs field of charge- 1 
and action S K . We will be interested in the small-K phase of the theory, when the 
unitary Higgs field fluctuates strongly and decouples at long distances (see^ for a 
recent study). We have included both Majorana (yh) and Dirac (y) gauge invariant 
Yukawa terms in the mirror action, in order to explicitly break all global symmetries 
of the mirror fermions (except of the U(l) to be gauged, as explained in Section 



3.3) 



When the unitary Higgs field acts essentially as a random variable (modulo cor- 
relations induced by small nonzero k and by fermion backreaction), the similarity of 



s We use ±, which is most commonly used in 2d (instead of L/R) to denote chiral fermion compo- 
nents. 
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the Yukawa interactions in ( 38 1 with the multi fermi interactions of ( 30 ) is in that 
integrating out the essentially random Higgs field generates all multi-fermion inter- 
actions consistent with the symmetries. Based on experience with strong- Yukawa 
expansions in theories with naive or Wilson fermions, it is expected that in the 
large-y, Gxed-h limit, there is a symmetric phase. 



The lattice action (3.2 1 completely defines the theory via a path integral over 
the charged and neutral fermion fields, the unitary higgs field, as well as the gauge 
fields. We will not perform the integral over the lattice gauge fields (but will study 
the variation of the partition function with respect to the small changes of the 
gauge background). Using this model, we hope to address some of the questions 



listed in Section 3.5.1 for example, while (1), (2), (3), and (4) will be addressed, 
it is clear that we can not answer question (5) about decoupling an anomaly-free 
representation, since the mirror fermions in the "1-0" model are in an anomalous 
representation. Thus, (5) — arguably one of the most important questions — is left 
for (near-) future study. 

4.1.2. Splitting the vectorlike partition function. The "splitting theorem." 

The partition function of the "1-0" model — as well as that of any vectorlike theory — 
can be split into "light" and "mirror" parts in any gauge background. To split 
the partition function, one uses the definite-chirality eigenvectors of 75 and the 
projectors P± on the corresponding spaces: 



■j 5 Ui = -m , %uii = Wi , (40) 

E ' ^+ = E ' 



a. = Yu iU \ , p+ = y Wr w\ = x-P- , (4i) 



where we treat u, w as columns and u\ as rows. We also use the eigenvectors of 
75 (the latter are independent of the gauge background) and the projectors P±: 

IbVi = Vi , j 5 U = -U , (42) 
P+ = J2 Vivj , P-=Y = 1 - P + ■ ( 43 ) 



Using (40 1, (42 1, a general Dirac field can be decomposed into chiral components 
af,a^ via the 75 eigenvectors, while the conjugate spinor field ^f x is decomposed 
into chiral components ~ct\ 7 a~ using the 75 eigenvectors, as follows: 



We now apply the split (44) to the fields ij) and \ of the "1-0" model, and note 
that only the charged eigenvectors (of both light and mirror fields) depend on the 
gauge background. The expansions ( [44) ) of the "mirror" fields are explicitly given 
below: 

x+ = ]T fa* , x + = ]T f + u\ [o] , V- = E «-*< > ^- = E ^- w l W ■ ( 45 ) 
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Clearly, expansions similar to (|45j) hold for the "light" fields as well: 



After substitution of (45 46), the partition function of the model (38) splits as 
follows: 

Z[A;y,h] = Z lighi \A\ x — I— x Z mirror [A;y,h] . (47) 
J[A\ 

Here Zu g ht[A] = det ||(uJ[A] • D[A] ■ Vj)\\x (similar determinant for the neutral light 
spectator X—) is the light sector partition function. The jacobian J is a product of 



jacobians for the charged and neutral sectors; see Section 2.2. of^and Section 4.1.5 
of this paper for details. The mirror partition function is given, more explicitly, by 
an integral over the charged mirrors (a_,5_), neutral mirrors (/3+, /3 + ), and unitary 
scalar field: 

Z mirror [A;y,h} = J d 2 a- d 2 /3 + dcj> e ~ s '™ , (48) 



where the mirror action from ( 38 ) is expressed in terms of the integration variables 



a_ , /?-)_ and the eigenvectors via ( 45 ) and d<j) denotes a path integral over the phases 



of 4>. The mirror-fermion integral is thus a determinant which includes the kinetic 



term and Yukawa terms from ( 38 1 and the mirror partition function is the average 



of the determinant over the random (in the disordered k — > phase) unitary field 



When A ^ 0, the mirror partition function Z mirror |48|) depends on the gauge 
background through the operators entering S m i rror (the Neuberger-Dirac operator 



and the associated projectors that appear in (38)) as well as through the gauge 



background dependence of the eigenvectors of 75 used to split the partition function 



{wi[A}, see @5J|). 

The dependence of general chiral partition functions, with an arbitrary chiral 
action, (e.g., our Z m i rror ) on the gauge background was studied irP^l There, an im- 
portant technical result was derived: that under an arbitrary variation of the gauge 
background, the variations of the mirror partition function due to the variations of 
the eigenvectors (here: iOj[A|) and the operators entering the action (the Neuberger- 
Dirac operator D[A] and the projectors P± depending on it, collectively denoted by 
O below) factorize. Explicitly, the "splitting theorem" states that for an arbitrary 
variation S of the gauge background, defined by the first equality below |^] the chiral 
partition function changes as follows: 

Z mirror\A + <5 A] _ \^ r x f \ , / ^ 



6\ogZ mrror [A] = log [ " ±p = E(<K ! ■ wA + (~ 50) , (49) 

^mirror \A\ ■ \OkJ 



'A more general chiral partition function depends on more than one set of A-dependent chiral 
eigenvectors, and a sum over their variations will appear on the r.h.s. of | |49[ |. We stress again 
that l |49| is valid for arbitrary chiral mirror actions (naturally, |4IJ[ ) gives the known result for the 
bilinear chiral actions appearing in (|26|). 
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where "(•)" denotes an expectation value calculated with the partition function 
Zmirror- The splitting theorerrP^ is valid for general chiral partition functions (e.g., 
the action does not have to be bilinear in the fermion fields) and is an important 
result, because: 

(f ) It isolates the anomalies from the details of the interactions of the chiral the- 
ory and manifestly realizes on the lattice the idea that anomalies are deter- 
mined only by the representation of the fields and not by the details of the 
Lagrangianl^l 

(2) The splitting theorem is indispensible^ in the calculation of mirror chiral 
current-current correlators in a perturbative expansion in the gauge field. 



In regard to (1) above, as already noted in Section 2.5 anomalies in chiral lattice 



gauge theories manifest themselves as topological obstructions to smoothly defining 



the basis vectors (or more precisely, the "measure current," J2i($ w l ' w i) m (49)) 
globally over the gauge field configuration spaced This implies that the splitting 
of the partition function in the "1-0" model into "light" and "mirror" sectors is 
singular,^ since every sector is anomalous on its own (but the full vectorlike theory 
partition function is not singular as a function of the gauge background A). However, 
because the singularity is topological, much like the location of the Dirac string of a 
monopole, it can be moved around the gauge field configuration space. In particular, 
a perturbative expansion around A = is insensitive to the singularity. Thus, the 
topological obstruction to making the splitting of the vectorlike partition function 
into "light" and "mirror" smooth over the entire gauge configuration space has 
no bearing on studies of the mirror spectrum at vanishing gauge background or 
infinitcsimally close to it. 



4.1.3. The phase structure of the "1-0" model at strong mirror Yukawa 
coupling 



The split partition function ( 48 ) allows the phase structure of the mirror sector of the 
"1-0" model to be studied^. The theory was considered in zero gauge background 
and in the limit of strong Yukawa coupling (y oo), where the mirror- fermion 
kinetic terms were dropped. As a function of h (the ratio of Majorana to Dirac 
Yukawa couplings), a nontrivial phase structure was found. For h — > and h > 1 the 
theory is in a strong-coupling symmetric phase with only short-range correlations of 
the /7(l)-charged order parameters. On the other hand, for h ~ 0.7, evidence for a 
transition to a "broken" (more precisely, algebraically-ordered phase) was seen. The 
transition appears to be in the Bcrezinskii-Kosterlitz-Thouless universality class, as 
evidenced by the change in the density of vortices^. 

On Figure^H[3j we show the behavior of the scalar field susceptibility as a func- 
tion of h. The susceptibility, which can also be thought of as the zero-momentum 
propagator of <f> x (or the square of the correlation length), is defined in the usual 
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Fig. 3. Scalar-field susceptibility l |50| of the mirror sector of the "1-0" model for k = 0.1, at 
•y — > oo. The dashed line indicates the susceptibilities for the pure XY model with the same k 
(undistinguishable, within errors for N = 4, 8, 16). Large errors at h = 0.7 and 0.8 are due to the 
sign problem at h < 1. 



way: 



N 2 \ 



E 



(50) 



where the brackets indicate an expectation value computed with the partition func- 



tion (48), on a square N x N lattice. As the figure shows, the correlation length is 



of order the lattice spacing for h > 1. The study of (50) as well as of other suscep- 
tibilities (composed from fermion bilinears) in refP^ show that the theory is in the 
strong-coupling symmetric phase for h > 1 and small k (at large k, on the other 
hand, the theory is in the "broken" phase). 



4.1.4. The polarization operator and transversality 

The spectrum of theories with GW fermions can be difficult to determine, particu- 
larly when the theory is in a strong-coupling phase. Numerical simulations are often 
restricted to states created by particular operators and do not lead to an exhaus- 
tive search for all possible light degrees of freedom that may exist. For example, 
irP3 a search for massless mirror states charged under the U(l) (would-be-gauge) 
symmetry, created by several classes of operators, yielded a negative result. This 
created a puzzle, since in the absence of such states, the long distance physics of the 
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"1-0" model would then be described by an anomalous unbroken gauge theory — a 
situation argued to be impossible without a violation of some key principle (e.g., 
unitarity). 

Clearly, in studies of mirror-fermion decoupling, we are particularly interested 
in states charged under the gauge symmetrjj^] (here: U{1)). The photon vacuum 
polarization operator, defined as: 

, , 5 2 \nZ\A\ 

n^(x, y )- 



can serve as a probe for all possible light degrees of freedom that are coupled to 
the gauge field. Below, we first study properties of the polarization operator of any 
theory with a gauge invariant partition function, i.e: 



with S LJ A fJ/ (x) 



kxZ[A + 8u>A] = lnZ[A] , 
-V M a; x , and V^w x = u> x +n ~ u x . This implies that: 
SlnZ[A] 



E 

//,.r 



5A„{x) 



V M w(a:) = 



Takin S of (53) gives: 



5\i\Z[A] 



(52) 
(53) 

(54) 



which, by expanding in A^ around A^ — 0, implies transversality of all n-point 
functions: 



Ev; 



S n Z[A] 

' 6A li (x)SA Ml (xi) . ..SA^^Xn-i) 



= 



(55) 



A=0 



For future use, note that, more generally, if O and F are any functions of the 
gauge field A, not necessarily gauge invariant, under the gauge transformation S u , 



A^A~ V^w, 



J w O(A)=^^%V / ,^(,i. 



SAJx) 



Therefore, the following identity: 



Ev; 



F ■ 



SO 



dAAx) J Suj(x) 



F-5^0 



(56) 



(57) 



is always true assuming only that O and F are smooth functions of the gauge field. 

The transversality condition (55) is only valid when it is applied to the full 
partition function of a theory where the gauge symmetry is respected, such as the 



u For anomaly matching arguments for global symmetries that will not be gauged, as in Section 
[5] we imagine first weakly gauging the relevant global chiral symmetry, finding the relevant polar- 
ization operator and then turning off the gauge field, as in 1511. 
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"1-0" model or any other vectorlike theory (we have assumed that Z[A] is a smooth 
function of A at least in the vicinity of A = 0; this is always a valid assumption 
in vectorlike theories provided that the action is smooth with respect to the gauge 
background since the fermion measure is well defined). 



4.1.5. Further properties and exact relations obeyed by the mirror IT MZ , 

Here, we present what will turn out to be (see the next Section) a lengthy derivation 
of the anomaly matching condition on the lattice for the simple "1-0" model. In the 
end of this Section, we also list some exact (i.e., independent of couplings) properties 
of the mirror polarization operator. 

The anomaly matching condition follows from the transversality condition ( 55 ) 
of the vectorlike theory and the possibility to (locally) smoothly split the partition 
function into light and mirror parts. Thus, as in (147]) , we have: 



In Z\A] = In Zugut [A] - In J[A) + In 

^mirror 

[A] , (58) 



where the Jacobian of the change of variables (44 1 is giverJ^, in the charged sector 
of the "1-0" model (ip±), by: 

J[A] = det||^(x) t u i (aj) t || det\\vi(x)t 3 (x)\\ . (59) 

Note that ||i;j(a:)£j(x)|| is a 2N 2 x 2N 2 dimensional matrix, with x indexing rows 
and i, j-columns and that there is a similar, but ^-independent, Jacobian arising 
in the neutral sector (x±)- From (581, it follows that the polarization operator (51 ) 
also splits into "light" and "mirror" parts: 

IV (.t, y) = H«f *(x, y) + II™ (as, y) . (60) 

We include the Jacobian contribution into the "light" sector H l ^ ht as a convention. 

Now, under an arbitrary infinitesimal change S v of the gauge background, the 
"light" partition function Zu g ht[A] = det(u][A] ■ D[A] ■ Vj) transforms (see^ for a 
detailed derivation) as: 

5 v ]nZ Hght [A] = TiiP+D-H.D) + J2( S v u ] ' u j) = ^{P+D^^D) + %[A],(6l) 

j 

where the currents j^ ,u '[j4] are defined as: 

j%[A] = ]T(Vj • u 3 ) , j™ [A] = • Wj ) . (62) 

j j 
In terms of the currents (62 1, the variation of the Jacobian isP^l 

6 7l lnJ[A}=j™[A} + MA}. (63) 



By combining ( |61[ ) and (63), we find that the "light" plus Jacobian contribution to 
the change of ( |58| under a gauge transformation is: 

^ ln ^T = ~ J " [A] - 1 2Z2»^**W , (64) 
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where j u denotes the measure current (62), restricted to a gauge variation of the 
background. Now we take of i64n and use the identity ( 57 ) to find: 



En 



SHZ Hght [A]J- 



w i) - ^ Xv l5xx[A] 



= -£v; 



■h 



'[A] - -t^W 



A] 



(65) 



We have introduced the shorthand notation that we will use freely hereafter: 

S 



5A^(x)> 



(66) 



where we suppressed the space-time index and understand that it is included in a 
single symbol /i, as long as no ambiguities arise. Expanding (65) around A^ = to 
linear order leads to: 



£vpi^f'(^) = -£v* ^j»L4] -y v tr% xx [A] , (67) 

showing that the "light" polarization operator is not transversal. There are two 
contributions to its divergence: the first term, proportional to the derivative of the 
"measure current," exactly cancels with the identical contribution of the "mirror" 
sector, as will become clearer later. 

The second term on the r.h.s. of (67), proportional to the derivative of the 
topological lattice field tr% xx (that this is a topological lattice field follows from 
cqn. (23), which shows that tT-js xx is an integer and thus can not change under 
a small variation of the gauge background) represents the anomaly of the "light" 
fermions. To make contact with the anomaly in the continuum, we note that the 
topological lattice field can be expressed as: 

1 



tr 75a 



Vlh»[A] 



(68) 



where = V ^A L ,(x) — V u A tl {x) is the field strength A^ and h^[A] is a gauge 
invariant local current, and ei2 = 1. An explicit form of /i M [A] can be obtained with 
some work. 

By the local smoothness of the "light" -"mirror" split, the fact that the full 
partition function of 1-0 model is gauge invariant implies the divergence of the 
"light" and "mirror" polarization operator cancel exactly. Hence, from (65), we 
conclude that 



E K* u ™ rror (x, 2/) = £V* + 2 Mr ^ ' (69) 

Since the terms that depend on the measure current in the "light" and "mirror" 
polarization operators cancel, they can be considered a lattice artifact due to the 
choice of the fermion measure. The last term in the above equation is more impor- 
tant since it can be expressed in terms of some correlation functions in the mirror 
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sector and therefore represents physically meaningful quantities. The fact that its 
divergence should agree exactly with that of the "light" sector appears to suggest 
that there always exist a certain light degree of freedom, even in the mirror sector, 
in order to contribute the appropriate anomalous divergence of n™ rror . 

We note that the conclusion that the mirror polarization operator is non- 
transverse, and that its divergence is proportional to the variation of the topological 
lattice field, depends only on the fermion content in the "mirror" sector and not 
on details on the mirror-theory dynamics, such as the strength of its couplings. We 



will argue in Section 4.2 that eqn. (69) is the lattice version of the t' Hooft anomaly 



matching condition in the GW formalism in 2d. 



Equations similar to (67 69 ) also hold for the non-transverse higher derivatives 
of the light and mirror partition functions. However, in 2d, to study the interplay 
between the anomaly and the light degrees of freedom it is sufficient to consider 
only the polarization operator in a trivial gauge background. In 4d, the trivial 
gauge background analysis of this Section would have to be extended to the three- 
point function in order to capture the effect of the anomaly. The mirror polarization 
operator would still be useful to detect massless mirror states in 4d. 

To end this Section, we now list some exact properties of polarization operators 
that hold for general chiral theories, in particular for our mirror theory; these are 
derived in Appendix A of refill All equations below refer to the polarization operator 
in x-space and not to their Fourier transforms; recall that we absorb the space-time 



indices into /i, v (see (66)). The definition of n™ rror = 8 a 8 v \ogZ mirror , see (60 
58), implies that n™ rror is symmetric due to local smoothness of Z m i rror . 

For the discussion here, the most important property is that the mirror polar- 
ization operator n™ rror ' can be decomposed into a part that is a measure current 
derivative and the rest: 

yrmirror r -w < ytrnirror I {7C\\ 

As shown irPS n™ rror ' is always a total derivative: 

j^mirror I ^ jjmirror I (71) 

and, in addition, Y\™- rror ' j s exactly gauge invariant. Therefore, proceeding as in 



the derivation of ( 65 1 , we find 



while with respect to the first index, we have 



y.jjmtrror ' = Q , (72) 



VJD? W ' = M. (73) 



Now, the total mirror n™ rror is symmetric, but 6 v j™ and n™ rror ' are separately 
not, but obey: 

ijrmirror I yfmirror l\ c -w , r -w 77 (7A\ 
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where JF^ is the curvature of the measure current, which is a known local functional 
of the gauge field whose divergence V* gives half the anomaly. These results 
imply that the symmetric and antisymmetric parts of n™ rror ' each contribute 
half of the anomalous divergence ( 73 ) . Since all properties listed in this Section 



hold independently of the mirror action, in particular of the strength of the mirror 



couplings, the verification of ( 72 ) , ( 73 ) , and ( 74 ) , in a numerical simulation at strong 



mirror couplings provides an important check on its consistency. 

Finally, we stress that in our numerical simulations, we calculate n™ rror ', as 
the measure-current part of n™ rror , see (70), is exactly the opposite that of the 
light theory. Obtaining an expression for n™ rror ' in terms of mirror-theory correla- 
tion functions, which can be computed using Monte Carlo methods, is a somewhat 
arduous task which is accomplished by repeated use of the splitting theorem (see 



Section 4.1.2); the calculation is explained in detail i: 



4.2. Anomaly matching and its possible solutions 

As already explained, in what follows, we focus on the n™ rror ' contribution to the 
mirror theory polarization operator. In fact, we will focus only on its symmetric part, 
since its antisymmetric part, as we discussed in Section |4.1.5| is of little physical 
interest — as it is independent on details of the interactions and depends only on the 



field content, see eqn. (74), it contains no information about the mirror spectrum. 



Here, we will study both the real and imaginary parts of the polarization op- 
erator of the mirror sector. We will show that its imaginary part has a nonlocal 
contribution that gives rise to the anomaly. The real part of the polarization oper- 
ator is a universal probe of the number and nature of the massless charged degrees 
of freedom and thus gives information of the spectrum. 

Let us first explain the nonlocality of the imaginary part of the polarization op- 
erator. For brevity, below we omit the superscripts (prime and mirror) and simply 

We work in momentum space and denote Fourier transfor- 



refer to U™ rror ' as ±V 



mation of n Mt/ as IL^. Let us first explore the consequence of equation (69 1, which, 



together with ( 70 ) and ( 68 1 , implies that for small momentum q: 

1 



2tt 



6„ x q A + 0(q z ) . 



(75) 



The rhs of (75 1 is of course local. One wonders if the usual argument in the contin- 



uum that it can not be the divergence of a local expression applies on the lattice 
too. A quick argument showing that it does is as follows. Let c = — i/(27r). Ignoring 
the real and transversal part of 11^^, we can rewrite ( 75 ) as the set of two equations: 



n 22 = -(c + iL7 2 ) 



<12 



(76) 



ftn = (c 



2l) qi 



(c-nr 2 )-- 
q\ 



where we used II 12 = II 21 . Locality of II n and II 22 would require that c — II 12 = 
Aqi + c + II 12 = Bq 2 + where A and B are arbitrary constants and dots 
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denote higher powers of momenta. This leads to: 

n 12 = -c + Bq 2 +0(q 2 ) , (77) 
n 12 = c-A qi + 0{q 2 ) , 

conditions, which are clearly incompatible. 



On the other hand, a nonlocal solution of the anomaly conditions (76 1 is given 
by IIi2 = c(g| — <Zi)/(<7i + <?!)• This is, indeed, the form of the non-transverse part 
of the polarization operator for an anomalous theory in the continuum, where the 
general solution of in 2d is given by: 

iVto = C(,») {tf - r) + c^ f'^ir^ ' , (78) 

where C\ and C2 are two arbitrary functions. 

We now enumerate the possible forms for C\ and Ci- As mentioned above, in 
Euclidean space the anomaly appears in the imaginary part of the polarization oper- 
ator only. It is well-known in the continuum that in a unitary Lorentz (Euclidean) 



invariant theory the zero-momentum singularity in the solution of ( 75 ) is due to 
either a massless Goldstone boson or a massless fermiorEH The massless scalar 
or fermion will, of course, also give a nonlocal, but divergence free, contribution 
to the real part of the polarization operator. If the GW formalism is to lead to 
any reasonable continuum limit, we would expect this conclusion to remain true on 
lattice. 

In order to see what nonlocal contributions to the real part of the polarization 
operator to expect, consider first a "Green-Schwarz" scalar theory in 2d Euclidean 
space: 

Z GS [A] = Jvri e / d2 *(-§ (d^-A,f +l £F 12 ) ( 79 ) 
Explicit computation leads to: 

The polarization operator is not transverse and its divergence is ig M ITg.g| A _ (q) 



}-£ up q p , identical to (75). The real part of (80) contains obviously the contribution 



2tt 

of a massless scalar; notice the shift of the coefficient k — > n + ~^t^ due to the 
anomalous Green-Schwarz term. If our mirror theory has a massless Green-Schwarz 
scalar, it would have to manifest itself by contributing to both the real and imagi- 



nary parts of the mirror polarization operator in agreement with ( 80 1 in the small 
momentum limit. 

Similarly, a massless mirror charged chiral fermion would also contribute to both 
the real and imaginary parts of II^„: 



fir I M-if^^rV-^^ (si) 
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where the real part of n M „ (q) of the chiral fermion is equal to one- half that of the 
Dirac fermion in the Schwinger model. Eqn. (81) is the small-momentum limit of 
the contribution of a free chiral GW fermion to the basis-vector independent part 
of the polarization operator, i.e., exactly what a massless mirror at y — would 
contribute to the full polarization operator of the vectorlike theory; see® 5 . Note 
that, as opposed to the Green-Schwarz scalar (80), the coefficients of the real and 
imaginary parts of the massless chiral fermion H^v arc the same|^] 

If long-distance unitarity was violated, one could imagine that the imaginary 
part of the polarization operator had a -j pole, providing the correct anomalous 
divergence, but the real part did not — as if there was no particle in the spectrum 
responsible to the anomaly. Our Monte Carlo simulations for the "mirror" sector did 
not find any evidence of such kind. Instead, it was found that unitarity is respected 
anywhere in parameter space, so that any chiral lattice theory with GW fermions 
(e.g., one with strong Yukawa interactions) appears always to lead to a physically 
reasonable field theory in the continuum limit. 



4.3. The symmetry arguments vs explicit lattice simulations 

Before we describe the results of the Monte Carlo simulation of the mirror sector of 
the "1-0" model, let us briefly explain the idea how this was done. 

In the continuum, the contribution to the Fourier transform of the real part of 
the polarization operator due to massless particles takes the form (78): 



IV(fc) = 2C 



k 2 



(82) 



where we have rescaled G\ = 2C for convenience and denoted the momentum by k. 
For a single massless degree of freedom, we have from (80): 

2Cgs scalar — ~ K (83) 



for a Green-Schwarz scalar, and for a single charged Weyl fermion from ( J8l| ) : 

2C ch . ferm . ~ ~ ~ -0.159 . (84) 

Z7T 



The constant C in (82) thus reflects the number and nature of massless degrees 
of freedom in the theory. To help identify the form of li^ and the coefficient C, 
we notice that the polarization operator ( 82 ) in 2d has a directional singularity as 



v We should note that in continuum 2d, the two realizations of anomaly matching and |80jl 
are equivalent by bosonization (seJ3H for the study of ano maly free 2d "Pythagorean" models as 
well as the more recent refP3). Our findings of Section |4.3| of massless fermions vs massless scalars 
in the strongly-coupled mirror sector of the "1-0" model should be really phrased as "in the lattice 
formulation that we use the results are naturally described as due to either a 'goldstone' or a chiral 
fermion." If we knew how to bosonize a lattice fermion with a Neuberger-Dirac operator, perhaps 
the equivalence could be extended to finite lattice spacing. 
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k — >• because: 

nH(^|^ = (7 ( 1 - cos2 ^ 

n 2 i(0)|^ o = -Csin2</>, (85) 

where <fr is the polar angle in the ki t 2 plane, so that approaches a different limit 
as k approaches from different directions. It is this angular singularity that can 
be easily identified in lattice simulations. In comparison, if only massive particles 
exist in the spectrum, the low momentum limit of II^,, would be proportional to 
(k^k 1 ' — k 2 g^ 1 ') I 'm 2 and the angular singularity is absent. 

Monte-Carlo simulations have been done for the mirror sector of the "1-0" model 
and the behavior the photon vacuum polarization operator in this theory has been 
investigated. The detailed results are presented irP^. It suffices to mention here 
that the system appears to adjust itself "automatically" in different regions of the 
parameter space, such that a single massless degree of freedom is always present, 
responsible for satisfying the anomaly equation. The angular singularity and the 



value of the discontinuity C ( 85 ) of the real part of Yl^ v that was found suggests 



that unitarity was respected. Furthermore, the value of C was found to agreep] well 



with either ( 83 1 or ( 84 1 , depending on the phase of the mirror theory, allowing us 
to interpret the result as due to either a "goldstone" boson or a chiral fermion. 
The agreement with the continuum values was impressive, given the rather small 8 2 
lattice used in the simulation (the complexity of the calculation of II M „ restricted 
us to rather small volumes). 

When the system is in the strong-coupling symmetric phase, where no massless 
scalar is expected, the result showed the consistency with the prediction of a single 
massless chiral GW fermion as shown in Figure |4j In this Figure, the behavior 
of II^j, with k approaching zero in three different directions is shown. The entire 
plot can be compared to the prediction of a single free GW chiral fermion, whose 
contribution to the divergence of can be analytically calculated. It was founcP^ 
that the two discontinuities agree very well. 

When the system is in the "broken" phase (at large k), where a single massless 
scalar should appear, the simulation also found the result to be consistent with the 
prediction given by a single Green-Schwarz scalar, as shown on Figure [5j 

Furthermore, when h — > in the symmetric phase, the discontinuity ( |85[ ) found 
in the Monte Carlo simulations was found to be approximately that appropriate 
for three massless chiral fermions. These are a massless charged chiral fermion and 
a massless charged vectorlike pair, whose appearance at h = can be analytically 



understood and will be explained in Section 4.4 

To summarize, the numerical evidence available so far suggests that strong 
Yukawa models formulated with GW fermions respect 't Hooft anomaly match- 

w The reader can consult ref.HU, or, alternatively, infer the small-fc discontinuities \85\ at <fi=0°, 
45°, 90° from the plots on Figs. |4| and [5| (replacing the index by 2) and then compare with the 
corresponding continuum values (|83| 1841 . 
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Fig. 4. The real parts of rioo and riio of the mirror for re = 0.1, h = 2, as a function of momentum 
approaching the origin in different directions. The value of the discontinuity indicates that when 
the system is in the strong-coupling symmetric phase, a single massless charged chiral fcrmion 
exists. 
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Fig. 5. The real parts of IIoo and Elio of the mirror for k = 5, h = 2 indicating the system is 
in a "broken" (algebraically ordered) phase where a single massless "goldstone" scalar appears, as 
explained i rM 



ing. Most importantly, in the strong coupling symmetric phase, where all mirror 
symmetries (but the one to be gauged) are explicitly broken, the minimal number 
of charged fermions required to saturate the 't Hooft anomaly of the to-be-gauged 
U(l) remains massless. One can view the results of this Section as a numerical 
"proof of 't Hooft anomaly matching and as evidence for the long-distance unitar- 
ity of GW-fcrmion Yukawa models. 

4.4. A different form of the split partition function; the relation to 
domain wall fermions 

In this Section, we will attempt to get a better analytic understanding of some of the 
numerical results of the previous Sections. To this end, we will develop a different 
representation of the 1-0 model partition function and comment on its relation 
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to domain-wall fermions. A representation like the one considered here, as well as 
the comments we make, hold for any GW-fermion model attempting to decouple 
mirror fermions from a vectorlike gauge theory by means of some non-gauge strong 
interactions. 



We begin by noting that the "1-0" model action (38) can be written in unpro- 



jected components (the superscripts P 1 ^ in the hatted projectors indicate whether 
the charge-1 or charge-0 Dirac operator has to be used): 

Skinetic = - (ipDilp) - (xA)X) 
S Y uka W a = U { (yPl<f>*P+x) + (XP-^P-^) } (86) 

+ y h{^ T (P-) T 4>j2P+x) - (xP-i2<i>*(Pl) T ^ T )} , 

where the measure now is the usual vectorlike theory measure in terms of ip, ip, \, X : 

dipdApdxdxd(f> e - s <— «<=- s ^" fc — " 5 « , (87) 



and S K is defined in ( 39 1 . We now perform the field redefinition: 



tp — > tp, tp — > tp 



1 



x -> x, x -> X 



1 



2-Dn 



The motivation for this redefinition can be traced back to the GW relation, which 
implies 75 = 75 93737 m other words, a 75 action on the original tp is trans- 



formed into the action of 75 on the new fields in ( 88 ) 



The price to pay for having an action of the lattice chiral symmetries generated 



by 75 j exactly as in the continuum, is the nonlocality of the redefinition (88) and 



thus of the resulting action (see (89) below). However, since we work perturbatively 
in the gauge field and at finite volume, we will imagine throughout this Section that 



the singularity at D — 2 of (88) is avoided by turning on background Wilson lines. 



The effect of the field redefinition on the partition function is: 



Si 



kinetic 



= -U> 



dtpdipdxdxdcp det(2 - D x ) det(2 - D ) e 
Di A (- D 



-SL 



2-D 1 



1> - X 



2 -D 



X 



(89) 



yh 



(iP T <p l2 P+x) 



XP-"f2<l>*tp 



xP--, 



Do 



Do 



12 <J 



D l 



(90) 



Obtaining the transformed Yukawa couplings in S' Yukawa requires repeated use of 



the GW relation and the equivalent relation shown in the paragraph after eqn. ( 88 1 



Needless to say, similar redefinitions hold in four dimensions and can be straight- 
forwardly performed if necessary. 
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Several comments, concerning the representation (89 90 1 of the mirror partition 
function, are now in order. We hope that these comments are useful to clarify the 
relation between different formulations of exact lattice chirality: 

(1) The singularity of the action at the position of the doublers, D = 2, ensures 
that they have infinite action and do not propagate at the classical level, as first 
proposed by RebbP^. The problem witfPSI, pointed out irJ 97 ! 98 !, of the doublers 
contributing as ghosts to the photon polarization operator at the quantum 
level is solved by the determinant prefactors, which exactly cancel the would-be 
ghost/doubler contributions. 



(2) Another comment concerns the relation of (89), with y = h = 0, to domain 
wall fermions. In the 2d case, these propagate on a finite interval in three di- 
mensions. The generating functional of Green's functions of the boundary chiral 
modes — in the case at hand, one charged and one neutral — can be obtained by 
integrating out the bulk fermions. This is technically possible since the fcrmion 
action is bilinear and there is no gauge field propagation in the extra dimen- 
sion (or any other non-uniformity except at the boundaries). When taking the 
chirally-symmetric limit of an infinite number of sites in the extra dimension, a 
massive bulk Pauli-Villars field, antiperiodic in the extra dimension, has to be 
included in order to obtain a finite determinant. The generating functional of 
Green's functions for the boundary chiral modes can be represented as a par- 



tition function with source terms. The result is exactly (89), with y = h = 0, 



and with source terms for ip and \ included. The determinant prefactors in 



( 89 1 arise as a combination of the determinants of the bulk fermions and the 
Pauli-Villars fields. The derivation of these results can be extracted from the 
work of refpSl. 

(3) It is less straightforward to relate the domain- wall fcrmion construction to mod- 
els where strong GW-fermion Yukawa or multifcrmion interactions are present. 
In this regard, recall the "pre-GW" proposal oP^, where strong multifermion 
interactions were added on one domain wall with the intent to decouple the 
mirrors from the Standard Model. We note also that the "warped domain wall" 
idea of refP^, whose 2d formulation was studied irP^, also uses strong interac- 
tions in a domain-wall set up. It would be interesting to more explicitly relate 
the finite-size extra dimensional version of these ideas (with strong "brane- 
localized" multifermion of Yukawa interactions) to the present (or similar) con- 
struction, since understanding the relation may yield both practical benefits 
and theoretical insight. 

(4) For vanishing Majorana coupling, h = 0, the Yukawa interaction in ([89]) is 



equivalent to the chirally invariant Yukawa coupling of refP^, see also^fil. To see 
this, use (ip^D^) = \k^D^) + \i^D Tj^Yji)) to replace the ip kinetic term in the 
action S' kin . Then note that -|(#D ^ V) = (?(2-D)0+^(?^)+^(^i50» 
where it is understood that the new charged field £ is integrated out from the 
action using its equation of motion, while the ^-determinant exactly cancels the 
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one in (89). Then, shift the integration variable ip — > rp + v2£. Next, perform 



exactly the same operations on x, introducing a neutral field 77, to finally obtain 
the action in the original form of Luscher, suitably adapted to the 2d case and 
to our normalization: 
1 



+ V20<I>*P + (X + \/2'7) + 77 (X + >/2rj)<l>P-(il> + V2£), 



27777 
V, 



(91) 



where the measure is the trivial one over ip, \, ?7, £• The Yukawa interaction of 
(3.2) is thus equivalent to that oP^l. We note that gauge and chiral invariant 



Majorana couplings were not considered irP2l and to the best of our knowledge 
were first constructed irP2 (see alscEO). 



For our purposes, the representation ( 89|90 1 of the "1-0" model partition func- 
tion is useful to emphasize the importance of Majorana couplings. It will allow us 
to understand the result mentioned at the end of Section [4.3| that at vanishing Ma- 
jorana couplings, the solution of the 't Hooft anomaly matching condition is not 
minimal and involves an extra vectorlike pair. The Majorana h ^ term gives mass 
to this extra massless vectorlike pair. 

We can split the partition function ( 89|90 ) exactly as we did in the basis of G W 
fcrmions. The representation of the split partition function that we give below is, in 



fact, equivalent to that in (47). The splitting of the partition function (89) into 
a "light" and "mirror" part is now done trivially using the 75-eigenvectors, which 
have no gauge-field dependence. We obtain, denoting now by ip± , x± the "normal" 
75-chirality components of the 2-component V'jX; an d using the fact that S' Yukawa 
only depends on the mirror components of ip,x : 

Z = Z+ " 1 



x Z. 



det_(2 



J ' 
A))det + (2 



Di 



di/j+dx- 



>4 



-V4 



where the mirror partition functiorj^ now is: 

Z_ = 

dct + (2 - D )dct_(2 - Di) ( dip^dx+d^e'' 



:x+)-s' Y 



(92) 

-S K 



Splitting the determinant prefactor into "light" and "mirror," as indicated in 
( 92|92 ), requires using the gauge-field dependent eigenvectors of 75: 

det+{2-D) = det + (l+757 5 ) = det||(ut(2-£%)|| = det||2(«tti)ll> ( 93 ) 
det_(2 -D) = det_(l +7575) = det||(w](2 - D)v 3 )\\ = det||2(wju,-)||, 
where the appropriate 75-eigenvectors are to be used for D° or D 1 . 



x Here and below, we renamed Z, 



light 



m 1 rror 



= z. 
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The gauge variation of the mirror partition function Z_ has now two contribu- 
tions: one from the variation of det_(2 — D 1 ) and one from the variation of the 
path integral over x+i 4>- The splitting theorem, applied to the chiral partition 
function defined by the path integral in ( 92 1 can be easily seen to imply that the 
gauge variation Z_ with the determinant factors left out vanishes when h = and 
y —> oo. Thus, the entire gauge variation of Z_ comes from the determinants: 

(51nZ_ = 5 a) lndet||(«;](2 — D 1 )v J )\\ = ^^(S u w^Wj) + tr P + —^——5 UJ (2 - D 1 ) 



(S u w] Wj ) + iTtuj (P{ - P+) = ^(Svwtwj) + ^Tr^, (94) 



3 3 



where, as usual, the measure current is cancelled by the variation of the Jacobian 
and light partition function. The gauge variation of Z_ is, naturally, the same as 
in ( 49 1 (the gauge variation of Z + can be obtained similarly to ( 94 ) and be seen to 
combine, together with the Jacobian to yield (64)). 

The determinant prefactor in Z_ contributes both to the real and imaginary 
parts of the mirror polarization operator. It is clear from (94 1 above that the imagi- 
nary part of the polarization operator due to the determinant is exactly as required 
by anomaly cancellation. Assuming unitarity at long distances, one would argue 
that Ren^„ should receive a contribution from at least one massless charged chiral 
fermion (plus, possibly, a number of massless states in anomaly-free representa- 
tions). The goal of our simulation was to precisely find out the real part of the 
mirror polarization operator. 

Let us now compare the numerical findings described in Section 4.3 with what 
can be inferred from the representation of the mirror partition function (92 92 1. 

The real part of the mirror theory 11"^ receives two contributions at y = oo. The 
first is due to the determinant prefactor. The contribution to Rell"^ of det_ (2 — 
D 1 ) can be seerP^to be exactly that of three massless propagating charged chiral 
fermions — the chiral components of the three 2d doubler modes. The determinant 
prefactor contribution to Ren~„ is, obviously, the same for any value of y, h. 

The only other contribution to Rell"^ at y = oo arises from the nonlocal coupling 
in the Majorana Yukawa term in S' Yukawa of (89). While it appears difficult to 
calculate analytically this contribution in the disordered-^ phase, our numerical 
results for the mirror polarization operator ReIT^„ show that there is one massless 
propagating chiral fermion. Thus two of the three massless modes contributed by the 
prefactor are cancelled by the nonlocal contribution to RelT^^ from the Majorana 
term in Syukawa', it appears that this cancellation is exact for all values of h > X. 

For vanishing Majorana coupling, h — 0, on the other hand, both numerical 
simulations (with the code of ^ used, we can only approach the h — > limit) and 
the representation ( 92 ) show that the real part of the mirror polarization operator 
at small momenta and infinite y is that of three charged massless chiral fermions. 
This is not too surprising as the Majorana coupling goes away in the limit and at 
infinite y the only dependence on the gauge field is in the determinant prefactor in 
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Z_ of (92 1. This indicates the crucial role of the Majorana-type couplings (recall 
that they were motivated by the need to break all mirror global symmetries) in 
facilitating the decoupling of the maximal possible number, as allowed by anomaly 
matching, of charged mirror degrees of freedom. 

Another lesson we learned is that probing the fermion spectrum with local 
fermion operators (including charged local fermion-scalar composites, as irP^) can 
miss massless degrees of freedom. The massless charged mirror fermions were not 
seen in that study, perhaps because they are not expressed in an obvious local way 
through the original variables. The long-distance gauge-boson polarization operator 
of the mirror theory is a universal probe of the charged mirror spectrum and should 
be the first quantity, along with susceptibilities probing chiral symmetry breaking, 
computed in any future studies of anomaly-free models. 



5. Outlook 

We have seen that when GW chiral fermions are used, the question of decoupling the 
mirror fermions from a vectorlike gauge theory is intimately intertwined with their 
contributions to anomalies. While this conclusion is not unexpected, these issues 
were usually not the central problem in earlier studies — largely because correlators 
of global chiral currents in the mirror theory could not be studied, as the exact 
lattice chiral symmetries were not known. This is where the GW formalism exhibits 
its advantages: when it is used, a precise formulation on the lattice of 't Hooft's 
anomaly matching condition can be found and verified. 

We have seen that models with mirror Yukawa interactions formulated via exact 
lattice chirality are "smart" enough that, consistent with anomaly matching, in 
the strong-coupling symmetric phase they lead to massless degrees of freedom as 
required by anomaly matching, rather than to a nonunitary theory. In the "1-0" 
model, the minimal number of chiral fermions needed to cancel the light fermion 
anomaly remain massless in the strong-coupling symmetric phase when all mirror 
global symmetries are explicitly broken. 

An obvious question is what these results imply for anomaly-free models. 

Take, for example, the "3-4-5" model, a two dimensional chiral U(l) gauge theory 
with 3_, 4_ and 5+ massless fermions, where the number indicates the U(l) charge 
and ± the chirality (we stick to 2d anomaly-free examples, as they are the logical 
next step; it is not realistic to expect the mirror dynamics of, e.g., the 4d SU(5) 
model of Section [3] to be studied numerically soon.) 

Consider first implementing the decoupling or mirrors from a vectorlike theory 
with three Dirac fermions of charges 3, 4, 5 by adding three uncharged mirror 
fermions and three unitary scalars, i.e. taking three copies of the "1-0" model and 
appropriately changing charges/chiralities (one could call the three copies the "3- 
0", "4-0", and "5-0" models). This is indicated on the lhs of Fig. [6] It is clear 
that such an implementation of the "3-4-5" model would lead to massless mirror 
states already in trivial gauge backgrounds, despite two facts that might suggest 
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0- 



3+\ / 4+ 



0'- 




/ 5- 
0+ 



3- 4- 5+ 
0+ 0'+ 0- 




3- 4- 5+ 0+ 



Fig. 6. Two possible ways to engineer mirror-fermion dynamics in the "345" model; the circles 
are meant to indicate that all interactions that break the mirror global symmetries between the 
enclosed fields are included. On the left: three copies of the "1-0" model leave three anomaly- 
free f/(l)'s whose 't Hooft anomalies require the presence of at least three charged masslcss chiral 
mirror fermions. On the right: the only unbroken global mirror symmetry is the to-be-gauged (7(1); 
since the mirror spectrum is anomaly free, 't Hooft anomalies do not require presence of charged 
massless states. Conjecture (testable!): the extra interactions breaking all global symmetries give 
mass to all mirror massless states in the strong-coupling symmetric phase — "just like" including 
the Majorana coupling in the "1-0" model lifted the massless vectorlike pair. 



otherwise: that the mirror and light spectra are anomaly- free and that, when the 
global U{1) appropriate to yield the "3-4-5" model is gauged, both the anomaly- 
free and anomalous global symmetries of such a lattice implementation are as in the 
target continuum theory. The easiest way to see that massless modes will result at 
zero gauge background is to weakly gauge the three separate global U{1) symmetries 
present in the "3-0" , "4-0" , and "5-0" models and demand consistency of the two- 
point current correlators, as we did for the "1-0" model above (the polarization 
operators for these U{l)s are essentially a set of correlation functions in the theory 
and we evaluate their divergences only at zero gauge-field background; thus all 
conclusions of Section [4] hold). 

The condition, already formulated in Section |3.3| but repeated here for com- 
pleteness, that the mirror interactions should obey in order to avoid massless states 
due to anomaly matching from an extra global symmetry is therefore to demand 
that when the gauge interactions are turned off, the mirror theory should have no 
global symmetries other than the global part of the gauge group — as elucidated in 
Section [373] With gauge interactions turned off, the implementation of the "3-4-5" 



Chiral lattice gauge theories via mirror-fermion decoupling: a mission (im)possible? 49 



model of the previous paragraph has, instead, two extra global f(l)s, which act si- 
multaneously on the light and mirror components and impose further conditions on 
the mirror spectrum that imply the existence of massless modes. This example leads 
us to conjecture that if the mirror interactions couple the 3+, 4 + , and 5_ mirrors, 
as schematically shown on the rhs of Fig. |6j by adding only one scalar and a single 
neutral 0_ mirror fermion and including the most general gauge-invariant couplings 
breaking all mirror global chiral symmetries, there wouldn't be any massless mirror 
states in the strong-coupling symmetric phase. 

At the moment, the strongest argument for the validity of our conjecture is that 
with all mirror global symmetries explicitly broken, there is no reason we know 
of, such as anomaly cancellation of any symmetry, for massless mirrors to exist at 
strong coupling (note that strong coupling is a must: since the gauge symmetry 
forbids mirror fermion mass terms, at weak mirror couplings and in an unbroken 
phase the mirror fermions will stay massless). The qualitative expectation at strong 
coupling is that the extra couplings on the rhs of Fig. [6] (compared to the lhs) that 
break all extra global symmetries will give mass to the massless mirror states that 
exist on the rhs (just like the Majorana coupling h of the "1-0" model gave mass to 
the vectorlike pair which remained massless at h = in the strong symmetric phase, 
as explained in Section 4.4). Needless to say, understanding how this mechanism 
works (or fails) in more detail would be very interesting. 

We end by noting that while achieving analytic understanding of the role of 
various couplings in the strong-coupling symmetric phase is quite desirable (but 
might be difficult), our conjecture about mirror-fermion decoupling can be tested 
via a numerical "experiment," similar to the one oP^, but this time for an anomaly- 
free model. Since no dynamical gauge fields are needed at this stage, the study of 
anomaly-free mirror dynamics is feasible given appropriate computer resources and 
the techniques already developed ir l 91 l 94 l 95 l and reviewed here. 
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